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I.  INTRODUCTION 


Armor-piercing,  fin-stabilized,  sabot-discarding  (APFSDS)  kinetic  energy  projectile 
configurations  consist  of  long  rod  cone-cylinder  bodies  with  a  set  of  teiil  fins  to  provide 
aerodynamic  stability.  An  edtemate  design  consists  of  a  flared  afterbody  as  the  sta¬ 
bilizing  element.  The  flared  projectile  design  has  been  investigated  for  training  round 
applications.^*^  Recently,  interest  has  been  generated  in  extending  the  flight  Mach  ntim- 
ber  for  APFSDS  projectiles  as  high  as  8.0.  Investigations  of  electro- magnetic  (EM)^  and 
electro-thermal  (ET)®  gun  propulsion  systems  as  well  as  ramaccelerator  (or  in-bore  ramjet) 
gun  systems®  have  shown  that  the  projectile  muzzle  velocity  could  reach  3  km/s.  These 
projectiles  generally  fly  low  altitude  trajectories  at  small  angles  of  yaw.  Under  these  con¬ 
ditions  the  effects  of  aerodynamic  heating  and  thermal  damage  to  the  projectile  become 
critical  design  issues.  The  need  for  computationally  rapid  projectile  design  codes  that 
can  accurately  predict  aerodynamic  drag  and  moment  coefficients  as  well  as  surface  heat 
transfer,  is  apparent.  These  codes  form  the  initizJ  stage  of  a  design  cycle  that  eJso  in¬ 
cludes  real-gas  (equilibrium-air  or  reacting-flow)  Navier-Stokes  fluid  dynamics  codes  that 
are  computationally  expensive. 

Preliminary  design  and  evaluation  of  proposed  projectile  configurations  require  rapid, 
reasonably  accurate  methods  to  predict  aerodynamic  characteristics.  Two  different  ap¬ 
proaches  have  predominantly  been  employed  to  calculate  aerodynamic  coefficients  using 
predictive  methods;  computational  fluid  dynamics  (CFD)  using  numerical  techniques  to 
solve  the  full  Navier-Stokes  equations  ^lnd  rapid  engineering  approximate  codes  using  semi- 
empirical  and  analytical  methods.  The  first  approach  is  more  generalized  and  can  accept 
a  wide  variety  of  freestream  conditions  and  body  geometries.  In  principle,  the  results  ob¬ 
tained  are  more  accurate  than  those  obtained  using  the  second  approach;  however,  CFD 
methods  require  considerable  computational  time  and  cost,  and  their  accuracy  depends 
on  selected  numerical  grid  sizes.  These  difficulties  have  restricted  CFD  applications  from 
engineering  problems  where  rapid  design  tools  are  required.  Instead,  CFD  serves  as  a 
followup  to  engineering  design  codes,  where  the  latter  is  used  to  rapidly  narrow  the  field 
of  viable  projectile  design  ceindidates  while  the  former  is  used  to  finalized  the  design  and 
assess  the  eierodynamic  impact  of  fine  geometric  changes. 

Engineering  design  codes  employ  approximate  linear  analyses  as  well  eis  experimental 
data.  Most  of  the  classical  techniques  have  been  proven  over  the  years  to  be  accurate  aero- 
dynsunic  predictive  methods.  However,  a  considerable  amoxmt  of  research  is  still  needed  to 
improve  and  extend  these  predictive  tools.  In  general,  the  areas  of  necessary  improvement 

^  Ctimint,  I.,  “Dn§  and  SiuHHtt  TndteBi  for  FUrt-SuHlixed  Projectile*,''  U.S.  Armt  BtllUtie  Reteoreh  Loiorotorf, 
Aierdeen  Prooiug  Gro%nd,  MD.,  report  in  preporution. 

^Celmini,  “Aerodpnomie  CAmneterutie*  of  Fin-  and  Flare-SiaHlixed  tSmm  XM910  Prototppe*,"  BRL-TR-MSt,  U.S. 
Armp  Balliitie  Reiearek  LaAoratorp,  Aberdeen  Prooinp  Gronnd,  MD.,  Dee.  1987. 

^Mermapen,  W.H.,  and  Yalamanekili,  R.J.,  “The  Effect  of  Perforation*  on  the  Balliotie*  of  a  Flare-StaAilized  Projectile," 
AIAA-8t-t0t4,  Proceedinp*  of  the  l3tA  AIAA  AtmoepAeric  FlipAt  Mechanic*  Conference,  Williamoharp  VA,  Aapaoi  18-90, 
1988. 

*  Gamer,  J.M.,  Zielinoki,  A.E.,  Jamioon,  K.A.,  “Deeipn  and  Teetinp  of  •  Ma**- Stabilized  Projectile  for  a  Small  Caliber 
Eleetromapnetie  Gan,"  BRL-MR-3744,  U.S.  Armp  Ballietic  ReeearcA  Laboratorp,  Aberdeen  Provinp  Gronnd,  MD,  April  1989. 

^ Penpellep,  R.,  and  Sweetman,  B.,  “FMC*  CAP  (Combnotion  Anpmented  Plaema)  Gnn,"  International  Defence  Reniew, 
Jnlp  1989,  pp.  813-815. 

* Hertzberp,  A.,  Bmekner,  A.P.,  and  Bopdanoff,  D.W.,  "Ram  Accelerator:  A  New  Chemical  Method  for  Acceleratinp  Pro¬ 
jectile*  to  UltrahipA  Velocitie*."  AIAA  Jonmal.  Vol.  tS,  No.  t,  Febmarp  1988,  pp.  195-903. 
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include  extended  scope  of  applicability,  ease  of  usage,  and  code  robustness. 

This  report  examines  an  engineering  design  code,  INTERACT,  written  at  the  U.S. 
Army  Ballistic  Research  Laboratory,^  that  employs  a  viscous-inviscid  interaction  pro¬ 
cedure  for  evaluation  of  axisymmetric  projectiles.  Separate  solution  methods  for  invisdd 
(Euler)  and  viscous  (boundary  layer)  flowfields  are  coupled  by  an  interaction  model  that  in¬ 
cludes  an  iterative  solution  procedure.  Inviscid  methods  include:  linearized  perturbation 
theory,  second-order  shock-expansion  theory,  modified-Newtonian  theory,  and  combina¬ 
tions  of  these  theories  for  certain  sections  of  the  projectile  geometry.  The  viscous  method 
employed  is  a  finite-difference,  second-order,  Keller-Box  scheme  for  standaurd  and  inverse 
solutions  of  the  boundary-layer  equations.  INTERACT  yields  body  surface  flow  profiles 
and  boundary  layer  profiles  in  less  than  one  minute  of  rxm  time  on  minicomputers.  These 
surface  profiles  represent  converged  solutions  to  both  the  inviscid  and  viscous  equations 
rather  than  an  inviscid  solution  updated  by  the  presence  of  the  boundary  layer.  The  ca¬ 
pability  of  computing  local  reverse  flow  regions  is  included.  The  procedure  is  formulated 
for  supersonic  and  hypersonic  Mach  numbers  including  both  laminar  and  turbulent  flow 
Reynolds  numbers.  Real-gas  and  high  temperature  effects  are  not  included  in  the  aerody¬ 
namic  analysis.  Approximate  methods  for  aerodynamic  heat  transfer  analysis,  assuming 
gas  at  chemical  equilibrium,  have  been  included.  The  heat  transfer  equations  use  body  sur¬ 
face  flow  distributions  and  compute  the  surface  heat  transfer  coefficient  and  local  Stanton 
number. 

Computations  of  surface  pressure  profiles  for  Mach  numbers  2  thru  5  and  8  are  com¬ 
pared  to  wind  tunnel  measurements  on  cone-cylinder-flare  projectiles.  Computations  of 
aerodynamic  stability  (center  of  pressure)  for  Mach  numbers  3  thru  6  are  compared  to 
wind  tunnel  measurements  on  a  cone- cylinder-flare  projectile.  Computations  for  surface 
heat  transfer  coefficient  are  compzired  to  results  obtained  firom  wind  tunnel  me^lsurements 
on  cone-cylinder-flare,  flat  plate,  and  blunt-cone  models  at  Mach  numbers  5  and  10. 


II.  BACKGROUND 


In  recent  years  engineering  design  codes  have  been  constructed  for  the  rapid  prediction 
of  aerodynamic  coefficients.  Two  classes  of  codes  have  emerged:  semi-empirical  methods 
2md  component  build-up  methods.  The  code  derived  in  this  report  addresses  a  third  class; 
viscous-inviscid  interaction  methods.  Some  of  the  available  codes  are  described  below  and 
summarized  in  Table  1. 

Semi-empirical  codes  are  formulated  from  a  database  of  experimental  and  test  data 
spanning  certain  ranges  of  Mach  and  Reynolds  numbers  and  for  certain  projectile  configu¬ 
rations.  A  specific  code  is  usually  limited  to  either  spin-stabilized,  small  L/D  projectiles  or 
non-spinning,  fin-  or  flare-stabilized  long-rod  configmrations.  Limitations  on  the  accuracy 
and/or  applicability  of  the  code  occur  whenever  the  flight  conditions  and/or  projectile 
geometry  exceed  the  range  of  the  database.  Exzunples  of  semi-empirical  codes  for  spin- 
stabilized  projectiles  consisting  of  a  conical  or  ogival  nose  section,  cylinder  mid-section,  and 

^ N%te*,  M.J.,  '‘Comp%ttUon*l  Atndfnamiet  for  Ariofmnttrie  Cone-Cflinder-Flart  Projeeiilet  i*  Sapersonie  Fli/kt," 
AIAA-iS-tgS7,  Proeeedingo  of  ikt  Ittk  AIAA  Atmoopktric  Fligkt  Mttkaniei  Conftnnet,  Snovmaoo  CO,  Appati  tS-tl,  198S. 


2 


boattail  or  flared  afterbody  include:  McDRAG,®  Donovan  zmd  Wood,®  Donovan,^®  and 
Morris.^^  For  fin-  and  flare-stabilized  long-rod  projectiles  the  codes  include:  Donovan, 
Nusca,  and  Wood,^®  Donovan  and  Grollman,*®  and  Donovan.^^  The  PRODAS  code  by 
Whyte  et.al.^^*^  is  formulated  to  predict  the  aerodynamics  of  both  spin-  and  fin-stabilized 
projectiles  and  includes  interior,  exterior,  and  terminal  ballistics  modules. 

Component  build-up  codes  are  formulated  using  classical  theories  for  the  prediction 
of  linear  aerodynamics,  and  semi-empirical  formiilae  for  special  cases  (e.g.  high  angle 
of  attack,  turbulent  flow).  In  these  schemes,  various  theories  or  semi-empirical  formulae 
are  applied  to  isolated  sections  of  the  projectile  configuration  and  then  combined  with 
predicted  interference  effects  to  yield  the  overall  aerodynamic  coefficients.  While  these 
codes  are  more  easily  extended  to  flight  conditions  outside  the  range  of  test  data,  often  each 
element  of  the  scheme  does  not  have  the  same  range  of  validity  as  other  elements.  Examples 
of  component  build-up  codes  include;  McCoy’s  combination  of  Vam  Dyke  hybrid  theory 
with  Van  Driest  boundary-layer  theory,'®  the  USAF  Missile  DATCOM  code,'®  and  the 
Naval  Surfaw:e  Wauffare  Center  Aerodynamics  Prediction  (NSWC-AP)  code.^®  Application 
of  the  DATCOM  and  NSWC-AP  codes  to  both  spin-stabilized  and  fin-stabilized  projectiles 
hais  been  evaluated  by  Vukelich  and  Jenkins,^'  Sun  amd  Cummings,®®  and  Mikhail.®® 

Viscou.i-inviscid  interaction  theory  is  essentially  am  extension  of  the  Euler /boundary- 
layer  theories.  Euler/boundaxy-layer  theory  uses  a  solution  procedure  for  the  inviscid 
(Euler)  flow  equations  to  determine  the  edge  velocity  boundary  condition  for  the  boundau-y- 
layer  solution.  Once  the  boundary  layer  has  been  computed,  the  inviscid  procedure  is  used 

*McCoTf,  R.L.,  “McDRAG  -  A  Computer  Program  for  Ettimating  the  Drag  Coeffteienit  of  Projeetileo,"  BRL-TR-tt9S, 
U.S.  Armp  Ballietie  Reteareh  Laboratory,  Aberdeen  Proving  Gronni,  MD,  Feb.  t9Sl. 

^Donovan,  W.F.,  and  Wood,  S.A.,  “Antomatie  Plotting  Rontinea  for  Eaiimaling  Aerodynamic  Properties  of  Spin  Siabilixed 
Projectiles  in  Flat  Fire  Trajectories  a(  BRL~MR-St04,  V.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving 

Groand,  MD,  Oct.  19SS. 

Donovan,  W.F.,  “Hypothetical  Zero  Yaw  Drag  Trajectory  of  Spinning  Projectiles  Between  A/=r5  and  M—IO,”  BRL-MR- 
3404,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Groand,  MD,  Nov.  1984. 

Morris,  M.A.,  “A  Compater  Program  to  Predict  the  Major  Aerodynamic  Coefficients  of  Conventional  Shell  and  Ballet 
Body  Shapes,"  RARDE  Memorandam  Report  9/81,  Jane  1981. 

Donovan,  W.F.,  Nasca,  M.J.,  and  Wood,  S.A.,  “Aatomatic  Plotting  Roatines  for  Estimating  Static  Aerodynamic  Proper¬ 
ties  of  Long  Rod  Finned  Projectiles  for  S<M<S,"  BRL-MR-SIZS,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving 
Groand,  MD,  Aag.  1981. 

Donovan,  W.F.,  and  Grollman,  B.B.,  “Procedare  for  Estimating  Zero  Yaw  Drag  Coefficient  for  Long  Rod  Projectiles  at 
Mach  Nambers  from  t  to  S,"  BRL-MR-S819,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Groand,  MD,  Mar. 
1978. 

Donovan,  W.F.,  “Aatomatic  Plotting  Aontine*  for  Estimating  Static  Aerodynamic  Properties  of  Flared  Projectiles  for 
t<.M<S,”  BRL-MR-3S7S,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Groand,  MD,  Mar.  1987. 

Whyte,  R.H.,  “SPIN-73,  An  Updated  Version  of  the  SPINNER  Compater  Program,"  U.S.  Army  Pieatinny  Arsenal  Tech¬ 
nical  Report  TR-4SSS,  November  1973. 

*•  Whyte,  R.H.,  Barnett,  J.R.,  and  Hathaway,  W.H.,  “FIN-81  A  Compater  Program  to  Predict  Fin  Stabilized  Projectile 
Aerodynamic  Coefficients  and  Stability  Parameters,"  General  Electric  i4rtnament  Systems  Dept.,  April  1981. 

^’’Fischer,  M.A.,  “PRODAS  User  Mtnaal  -  Veraion  3.1,"  General  Electric  Xrmament  Systems  Dept.,  Febraary  1989. 

McCoy,  R.L.,  "Estimation  of  the  Static  Aerodynamic  Characteristics  of  Ordnance  Projectiles  at  Sapersonic  Speeds,"  BRL 
Report  No.  188t,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Groand,  MD,  Nov.  1973. 

**  Vakelieh,  S.R.,  and  Jenkins,  J.E.,  “Missile  Dateom:  Aerodynamic  Prediction  of  Conventional  Missiles  Using  Component 
Baild-Up  Teeknigaes,"  AIAA-84-0388,  1984. 

^ Devan,  L.,  and  Mason,  L.A.,  “Aerodynamics  of  Tactical  Weapons  to  Mach  Namber  8  and  Angle  of  Attack  180":  Part  II, 
Compater  Program  and  Users  Gaide,"  NSWC  TR  81-3S8,  Naval  Sarfaee  Warfare  Center,  Sept.  1981. 

Vakelieh,  S.R.,  and  Jenkins,  J.E.,  “Evalaation  of  Component  Baildap  Methods  for  Missile  Aerodynamic  /’reliction,* 
AIAA  Joarnal  of  Spacecraft  and  Rockets,  Vol.  19,  No.  8,  Nov/Dee  198i,  pp.  4dl-4d8. 

San,  J.,  and  Cammings,  R.M.,  “Evalaation  of  Missile  Aerodynamic  Characteristics  Usin,  Rapid  Prediction  Teehnigaes," 
AIAA  Joamal  of  Spacecraft  and  Rockets,  Vol.  tl.  No.  9,  Nov/Dee  1984,  PP-  513-StO. 

Mikhail,  A.G.,  “Interactive  Inpat  for  Projectile  Fast  Design  Codes  Package,"  BRL-MR-3631,  U.S.  Army  Ballutie  Research 
Laboratory,  Aberdeen  Proving  Groand,  MD,  October  1987. 
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to  recompute  the  flow  over  the  original  body  combined  with  the  botmdary  layer  displace¬ 
ment  thickness.  The  procedure  is  uncoupled  and  does  not  recompute  the  botmdary  layer. 
A  version  of  Euler /boundary- layer  theory  has  been  applied  to  spin-stabilized  projectiles 
in  transonic  and  supersonic  flight  by  Sturek,  Nietubicz,  and  others.^^^  The  code  by 
Sturek  et.al.*^  computes  the  inviscid  (Euler)  flow  over  a  projectile  using  a  finite-difference 
MacCormack  technique.  The  effect  of  the  boundary  layer  on  the  flowfield  adjacent  to  the 
surface  of  the  body  is  accotmted  for  by  adding  the  local  displacement  thickness  to  the 
local  body  rsKiius  and  solving  the  inviscid  flowfield  with  an  adjusted  boundary  condition. 
The  present  code,  INTERACT,  includes  an  interaction  model  for  the  viscous  and  inviscid 
theories  allowing  the  independent  solutions  to  be  coupled.  Thus,  the  resulting  flowfield 
represents  a  solution  to  the  botmdary-layer  equations  and  the  Euler  equations  simulta¬ 
neously.  The  use  of  an  inverse  solution  to  the  botmdary-layer  equations  (i.e.  specified 
boimdary  layer  displacement  thickness  and  computed  edge  velocity)  is  a  necessary  element 
of  the  procedure.  Similar  techniques  have  been  successfully  used  for  aiifoils  in  low  speed 
flow.*^'^  The  application  of  this  method  to  supersonic  flow  is  outlined  by  Nusca^  and 
described  in  Section  III.3  of  this  report. 

Aerodynamic  heating  analysis  for  projectiles  flying  at  Mach  number  5  was  examined 
by  Sedney^*  in  1958.  This  analysis  assumed  that  the  projectile  geometry  can  be  repre¬ 
sented  by  two-dimensional  components  for  which  approximate  values  of  the  heat  transfer 
coefiicient  have  been  documented  by  Van  Driest.^^  This  aneilysis  yielded  the  surface  tem¬ 
perature  versus  time  of  flight  (or  range)  for  small  L/D,  spin-stabilized  projectiles.  Since 
then  multi-dimensional  analyses  and  rigorous  experimental  tests  have  yielded  closed  form 
formulae  for  the  calculation  of  heat  transfer  coefficients  on  general  axisymmetric  geome¬ 
tries.  These  formulae  relate  surface  heat  transfer  to  computed  (or  measured)  boundary 
layer  profiles  and  boundary  layer  edge  flow  properties  (i.e.  velocity,  pressure,  entropy). 
Some  examples  of  this  work  include  Rubesin  and  Johnson,^  Young  and  Janssen,^  Eckert,^® 
Fay  and  Kemp,^  Crabtree  et.zd.,^^  Fancett,^  DeJarnette  et.al.,^®  and  recently  Hoffmann 
et.ad.^°  In  pau-ticular,  Fancett^  aind  Hoffmann^®  have  assembled  engineering  heat  transfer 
formulae  and  applied  these  methods  to  axisymmetric  projectile  configurations.  The  work 
of  Famcett^  is  used  here  (Section  IV)  to  compute  surfaice  heat  transfer  coefficients  from 

Sturek,  W.B.,  Dwyer,  H.A.,  Kayeer,  L.D.,  Nietehiez,  C.J.,  Reklit,  R.P.,  and  Opalka,  K.O.,  “Comp%taiion»  of  Moynxe 
EBoeU  for  a  Yawed,  Spinning  Body  of  Revolntion,"  AIAA  Journal.  Vol.  IS,  No.  7,  Jnly  1978,  pp.  S87-S9t. 

**S<«rei.,  W.B.,  Mylin,  D.C.,  and  Bath,  C.C.,  “Computational  Parametric  Study  of  the  Aerodynamiea  of  Spinning  Slender 
Bodiea  at  Superionie  Speedi,"  AIAA-80-IS8S,  Proceeding!  of  the  7th  AIAA  Atmoipherie  Flight  Mechanic!  Conference,  Augmt 
tl-tS,  1980. 

^ Nietukicx,  C.J.,  Dankerg,  J.E.,  and  Inger,  G.R.,  “A  Theoretical  and  Experimental  Inueatigation  of  •  7Vam«a«ie  Projectile 
Flow  Field,"  BRL~MR~St9l,  U.S.  Army  Balliatie  Reoeareh  Lakoratory,  Akerdeen  Proving  Ground,  MD.,  July  198S. 

^^Pleteher,  R.H.,  “Prediction  of  Ineompreaikle  Turkulent  Separating  Flow,"  Jonmal  of  Fluid!  Engineering,  Vol.  100,  Dee. 
1978,  pp.  dtf-dSJ. 

^*Kwon,  O.,  and  Pleteher,  R.H.,  “Prediction  of  Ineompreaikle  Separated  Boundary  Layer!  Including  Viecoue-Invieeid  In¬ 
teraction,"  Tran!aetion!_ofjkejlSME,  Vol.  101,  Dee.  1979,  pp.  dff-dH- 

^  William!,  B.R.,  “The  Prediction  of  Separated  Flow  Uoing  A  Vioeouo-Inviaeid  Interaction  Method,"  Aeronantieal  Journal. 
May  1985,  pp.  185-197. 

®*/m,  B.J.,  “A  Coupled  Vioeid-Invioeid  Calculation  Method  for  the  Prediction  of  Umteady  Separated  Flow!  Over  An  Airfoil," 
AIAA-88-05S8,  Proceeding!  of  the  tSth  AIAA  Aeroepace  Science!  Meeting,  Reno  NV,  Jan.  11-14,  1988. 

Sedney,  R.,  “Aerodynamic  Heating  Prokleme  In  Shell  Deeign,"  BRL  Report  No.  1043,  U.S.  Army  Balliatie  Reaeareh 
Lakoratory,  Akerdeen  Proving  Ground,  MD,  March  1958. 

**  Vkm  Drieet,  E.R.,  “The  Proklem  of  Aerodynamic  Heating,"  Aeronautical  Engineering  Review,  Vol.  15,  No.  10,  Oetoker 
1958,  pp.  96-41. 


4 


laminar  and  turbulent  flow  profiles  as  determined  by  the  INTERACT  code. 


T^ble  1.  AERODYNAMIC  DESIGN  CODES  FOR  AXISYMMETRIC  BODIES 


Code 

Inviscid  Methods 

Boundary  Layer 
Methods 

Inviscid/BL 

Interaction 

Heat  IVansfer 

MCDRAG* 

Similarity  rules 

Flat  plate 
empirical 

None 

None 

Donovan,  Wood® 

Empirical  data 
curve  fits 

Flat  plate 
empirical 

None 

None 

Donovan  (flare) 

Empirical  data 
curve  fits 

Flat  plate 
empirical 

None 

None 

McCoy‘* 

Van  Dyke  hybrid, 
Linear  potential 

Van  Driest  flat 
plate  theory 

Single 

iteration 

None 

DATCOM^® 

Component  build-up 
Slender  body  theory 

Flat  plate 
empirical 

None 

None 

NSWC-AP2° 

Van  Dyke  hybrid, 
2nd-order  Shock/Exp, 
Newtonian  theory 

Van  Driest  flat 
plate  theory 

None 

None 

INTERACT 

Linearized  perturbation, 
2nd-order  Shock/Exp 
with  flare  corrections, 
Newtonian  theory 

2D/Axisymmetric 
2nd-order 
finite  difference 
solution 

Multiple 

iteration 

w/convergence 

criteria 

Algebraic, 

assumes 

chemical 

equilibrium 

(No  attempt  made  to  include  all  available  codes) 

M.W.,  tnd  Johnson,  H.A.,  "A  Critical  Reoieiu  of  Skin- Friction  and  Hcat-Tran$/tr  Sointiont  of  the  Laminar 
Bonndarf  Laftr  of  a  Flat  Plate, ’'Tranoactiono  of  the  ASMS,  Vol.  71,  ^o.  4,  Afaf  1949,  ff.  3SS-38S. 

**  Yonng,  G.B.,  and  Janoocn,  E.,  “The  ComprtieHU  Bonndarf  Lafcr,"  Jonmal  of  Aeronantical  Scienceo,  Vol.  19,  April 
1959,  fp.  999-935. 

Eckert,  E.R.,  "Engineering  Relatione  for  Friction  and  Heat  Transfer  to  Snrfaees  in  High  Veloeitf  Flov," 
Jonmal  of  Aeronantical  Sciences,  Vol.  99,  No.  9,  Mvfwl  1955,  pp.  585-587. 

**/’•»,  J.A.,  and  Hemp,  N.H.,  “Theorp  of  Stagnation-Point  Heat  Transfer  in  a  Partiallp  Ionised  Diatomic  Gas," 
AIAA  Jonmal.  Vol.  1,  No.  19,  Dec.  1953,  pp.  9741-9751. 

Craitree,  L.F.,  Dommett,  R.L.,  and  IVoodlep,  J.G.,  “Estimation  of  Heat  Transfer  to  Flat  Plates,  Cones  and  BInnt  Bodies," 
RAE  TR  No.  3537,  Jnlp  1955. 

^Fancett,  R.K.,  "Heat  Transfer  to  Slender  Projectiles  at  Mack  Nnmkers  Between  5  and  10,"  RARDE  Memomndnm  Report 
9/70,  1970. 

^*DeJamette,  F.R.,  Hamilton,  H.H.,  Weilmnenster,  K.J.,  and  Ckeatwood,  F.M.,  “A  Review  of  Some  Approrimate  Methods 
Used  in  Aerodpnamie  Heating  Analpses,"  AIAA  Jonmal  of  Thermophpsics,  Vol.  1,  No.  I,  Jan.  1987,  pp.  5-19. 

Hoffman,  K.A.,  Wilson,  D.E.,  and  Hamhnrger,  C.,  MeT«ficrm»4|rnsmie  Analpses  of  Projectiles  at  Hppersenie  Speeds," 
AIAA-9185,  Proceedings  of  the  7th  AIAA  Applied  Aeredpnamies  Conference,  Seattle  WA,  Jnlp  31  to  Angnst  9,  1989. 
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III.  INTERACTION  MODEL 


1.  Inviscid  Methods 


The  inviscid  flowfield  is  governed  by  the  Euler  equations  (i.e.  Navier-Stokes  equations 
with  zero  kinematic  viscosity  or  infinite  Reynolds  number)  which  describe  the  conservation 
of  mass,  momentum  and  energy  for  a  compressible  flow. 

^+v.(^i?)  =  o  (1) 


DV 
^  Dt 


=  -Vp 


(2) 

(3) 


These  equations  can  be  solved  by  finite'difference  techniques  requiring  extensive  com¬ 
putational  resources.  However,  finite-difference  methods  are  not  the  preferred  choice  for 
aerodynamic  design  codes  since,  for  general  geometries,  grid  generation  requires  skilled 
interaction  between  the  code  and  operator.  Previous  aerodynamic  design  codes  have  uti¬ 
lized  linearized  perturbation,  shock-expansion,  and  modified-Newtonian  theories  to  solve 
the  inviscid  flowfield  (see  Table  1).  These  methods  axe  suitable  choices  for  the  present 
code,  since  they  yield  accurate  solutions  when  applied  within  the  limits  of  each  theory,  do 
not  require  grid  generation,  and  run  for  less  than  one  minute  CPU  time  on  minicomputers. 
Linearized  perturbation  and  shock-expansion  theories  are  limited  to  pointed  slender  bod¬ 
ies,  while  the  limits  on  freestream  Mach  number  and  flow  deflection  angle  for  each  theory 
axe  quite  different.  Lineznized  perturbation  methods  axe  based  on  linearized  solutions  of 
Equations  1-3  written  in  velocity  potential  form.  Shock-expansion  theory  is  based  on  an 
eissumption  of  locally  two-dimensional  flow  with  a  correction  for  axisymmetric  effects.  For 
a  flared  afterbody,  a  special  correction  is  required.  Linearized  perturbation  theory  is  more 
generally  applicable  to  irregular  body  shapes,  but  imposes  greater  restrictions  on  the  body 
fineness  ratio  and  freestream  Mach  number  them  does  shock-expansion  theory.  Modified- 
Newtonian  theory  can  be  applied  on  the  forebody  section  of  projectiles  with  blunt  cone  or 
bltmt  ogive  geometries. 

Based  on  the  freestream  Mach  number,  the  INTERACT  code  employs  the  most  suit¬ 
able  method  or  combination  of  methods  for  a  given  projectile  geometry,  thus  accotmting 
for  most  cases  of  interest  to  the  projectile  designer.  These  methods  are  reviewed  in  the 
following  subsections. 


a.  Linearised  Perturbation  Theory 

Assiiming  that  the  flowfield  over  the  body  is  steady  and  irrotational,  then  ^  =  0  and 
V  X  F  =  0.  It  can  be  shown,^^  that  Equation  2  reduces  to 

dp  =  -pVdV  (4) 
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The  Irrotational  flow  assumption  is  exactly  true  for  the  flow  behind  an  attached  oblique 
shock  wave  on  a  wedge,  and  approximately  true  for  a  sharp-nose  slender  body  in  supersonic 
flow.  If  the  slender  body  is  moving  supersonically,  the  attendant  shock  wave  will  be  slightly 
curved,  and  hence,  strictly  speaking^  the  flowfield  will  be  slightly  rotational.  However,  it 
is  usually  practical  to  assume  V  x  V  0  for  such  cas^.  It  is  then  convenient  to  define  a 
velocity 'potential,  $,  such  that  V  =  where  V  x  V  =  V  x  V#  =  0,  with  #  a  scalar. 
Application  to  the  continuity  equation  (Equation  1)  and  to  Euler’s  equation  (Equation  4) 
for  isentropic  fiow  (i.e.  dp/ dp  =  a^),  yields  the  Velocity  Potential  Ekjuation: 


$2 


$2 


-|-(1-— -2—^ 
w  '  V-*-  .,2  ^  - 


— -  2-V-$ 


a* 


V* 


=  0  (5) 


where 

=  Of  ~  ^($1  +  the  local  speed  of  sound 

n  =  the  total  speed  of  sound,  constant  for  a  calorically  perfect  gas 


The  Velocity  Potential  Equation  (Equ.  5)  is  a  non-linear  partied  differentied  equation 
that  is  exact  within  the  framework  of  irrotational,  isentropic  flow.  It  applies  to  all  values  of 
the  Mach  number.  Equation  5  can  be  applied  to  a  slender  body  in  uniform  flow  where  the 
freestream  velocity  (Kw)  is  only  slightly  disturbed  (perturbed).  The  analysis  of  this  type 
of  flowfield  is  referred  to  as  small-perturbation  theory.  The  body  causes  a  disturbance  of 
the  uniform  flow  and  the  velocity  is  written  in^terms  of  the  perturbed  velocity  components, 
UjU,  tS,  in  the  x,y,z  directions  where  V$  =  V  =  (Voo  +  +  {^)j  +  (w)fc‘  The  perturbed 

velocity  potential,  is  defined  using: 

$(x,  y,  z)  =  VooX  +  4(z,  y,  z) 


Substituting  this  expression  into  Equation  5,  we  obtain  the  velocity  potential  equation 
for  the  perturbed  fiow: 


=  Ml 


(t  +  I)  .  (7-1) 


+ 


('”"c 


(^)) 


du 

dx 


-i-  Ml 


+  Ml 


(7  +  1)  (7-1)  ( ^ 

^  2  2  \ 

.  (7  +  1)  ,  (7-1)  ( “1+^ 


JJdy 

))i 


W 


(6) 


where 


**  Andtr$on,  J.D.,  Modern  Comprt$$Ule  Flew  with  Hittoricml  Pertpeeiiwe,  MeGnw-HiU,  198t. 
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a’  =  al,  -  ^(2uVoo  +  +  w^) 

Mfjo  —  V<x>/ Ooo 


Assuming  that  u,t5,and  w  are  small  compared  to  Voo,  then  u/V'oo,v/V’oo,t5/V’oo  <C 
1  and  (u/K»)*,(v/V'oo)^»(dj/V'oo)^  -C  1.  We  can  therefore  exclude  from  ^uation  6  the 
term  M^(7  +  l)(u/Voo)5fi/5i  in  favor  of  (1  —  M^)5u/5x  except  when  Moo  —  1  and  the 
terms  M^{‘y —  l)\u/Voo)dv/dy  and  M^{'f  —  l){tL/Voo)div/dz  in  favor  of  dv/dy  and  dwjdz, 
respectively,  for  Moo  <  5  (approximately).  We  then  obtain  the  Prandtl-Glauert  equation: 


dx^  ay2  '  dz^ 


=  0 


(7) 


or  in  axisymmetric  form  (iw  =  0), 


5r^  r 


(8) 


where,  0  =  -  1 

Equation  8  is  a  linear  equation  that  is  valid  for  inviscid,  irrotational,  isentropic  flow  with 
small  perturbations  to  the  freestream  flowfield.  In  addition,  flows  for  which  .8  <  Moo  ^  1*2 
and  Mao  >  5  (approximately)  have  been  excluded.  Von  Karman  and  Moore^^  showed  that 
Equation  8  can  be  solved  by  placing  a  source  distribution  along  the  axis  of  the  body  and 
using  the  method  of  superposition  of  elementary  solutions.  This  method  is  restricted  to 
axisymmetric  flow  (i.e.  the  body  is  at  zero  incidence  and  w  =  0).  The  source  strength  of 
each  element  is  located  at  i  =  ^  +  0r.  Integration  over  all  sources  for  the  entire 

length  of  the  body  (z  =  0  to  L)  yields  the  following  relation  for  the  perturbed  velocity 
potential: 

4(x,r)  =  /  fm - (9) 

where,  r  =  v/y*  +  2*;  the  radial  distjince  from  the  zixis  of  the  body. 

Equation  9  applies  to  a  pointed  body  where  /(^)  =  0  for  ^  <  0  (x  <  0).  The 
change  in  the  limits  of  integration  is  a  consequence  of  the  imaginary  contribution  to  $ 
\i  ^  >  X  —  0r  (i.e.  the  influence  of  sources  downstream  of  x  is  not  considered  for  ptirely 
supersonic  flow).  Using  the  transformation  ^  =  x~0r  cosh <7,  =  ^x  —  —  0^r^da  and 

applying  the  method  of  superposition  of  elementary  solutions  (i.e.  /'(O  =  Aj  =  constant 
for  i  <  (i,  dividing  the  body  axis  into  n  simple  parts  (t  =  1, . . .  ,n)  and  imposing 

separate  sources  in  these  parts),  the  velocity  components,  5,5,  are  then  expressed  as^^ 

“(®M  n)  ^cosh'*  -  cosh"^  (10) 

H.W.,  tnd  Rotko,  R.A.,  EUmtnt$  of  Wiltf,  NY,  1957  (ff.  ti7-iS0). 
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Combining  Equations  10  and  11  with  a  tangent  flow  boundary  condition  on  the  body 
surface,  drifdx;  =  v,7(K»  +  «»)  yields  a  system  of  three  linear  equations  for  the  unknowns 
u(x,,r,),  v(xi,ri),  and  Ai.  Since  Ai  is  a  function  of  Aj  for  j  =  l,...,i  —  1,  the  solution 
is  marched  from  Ai  =  0  at  i  =  0  to  x  =  x,-.  The  solution  must  be  stsjted  from  the 
nose  of  the  body.  Application  of  Equations  10  and  11  for  the  solution  of  the  inviscid  flow 
over  a  slender  body  has  been  termed  linearized  perturbation  theory.  The  mathematics  of 
Equations  10  and  11  dictate  the  following  condition  for  the  flow  deflection  caused  by  the 
body  nose  angle  (^„).  For  t  =  2: 


«2  =  A2(-cosh  ^[(X2  -  6)/(;^r2)]) 

V2  =  -Ai  (-i/(x2-6)^ 

Therefore,  (12  -  ^i)/(j3r2)  >  1  and  (x2  -  ^1)^  -  >  0  and  thus 

^  tan  ~  ^ 


(12) 


Since  the  solution  must  be  started  from  the  nose  of  the  body,  Equation  12  determines 
the  range  of  applicability  of  the  theory,  assuming  that  the  flow  remains  attached  to  the 
body.  Equation  12  is  equivalent  to  the  limit  S„  <  =  arcsin(  1/Moo)  (i-e.  that  the 

freestream  Mach  number  lines  reside  outside  the  body  surface).  Figure  1  displays  this 
limit  for  1  <  Moo  <  10. 


Other  body  surface  flow  properties,  Cp.  (local  pressure  coefficient),  Vi/V^o  (local  tan¬ 
gential  velocity),  and  Pi/Poo  (local  flow  density)  can  be  obtained  from 


Pi  _ 

Poo  P.oMi'VP 


(14) 

(15) 


where,  M„  and  p„  are  the  Mach  number  emd  pressure  behind  the  nose  shock  (see  Equ  18 
and  19). 


b.  Second-Order  Shock-Expansion  Theory 

Solutions  of  the  inviscid  flow  about  simple  geometries  can  be  obtained  easily  if  the 
flowfield  consists  entirely  of  shock  waves  and/or  expzmsion  waves.  For  generalized  (i.e. 
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first-order)  shock-expansion  theory,  the  flow  over  an  axisymmetric  body  at  zero  incidence 
to  the  fireestream  is  assumed  to  be  locally  two  dimensional.  Therefore,  the  continuity 
equation  (Equation  1)  is  only  approximately  satisfied.  Shock  waves  are  assumed  to  be 
oblique,  straight,  and  attached  to  the  body.  Body  surface  properties  of  the  fiow  (pressure, 
Mach  ntimber,  and  density)  can  be  obtained  by  application  of  the  oblique  shock  relations 
for  the  flow  behind  a  shock  wave  attached  to  the  nose  of  the  body,  and  the  Prandtl-Meyer 
function  for  compressions  and  expansions  caused  by  changes  in  body  slope  downstream  of 
the  nose.  Eggers  et.al."*^  documented  the  generalized  shock-expansion  theory  for  supersonic 
flow.  Second-order  shock-expansion  theory  was  proposed  by  Syvertson  and  Dennis^  to 
correct  the  first-order  resvilts  for  non-isentropic  and  axisymmetry  effects. 

In  general,  an  oblique  shock  wave  is  three-dimensional  and  curved.  However,  it  is 
usually  sufficient  to  consider  the  flow  at  an  arbitraiy  point  to  be  locally  two-dimensional 
and  the  shock  to  be  straight.  It  is  then  possible  to  derive  the  oblique  shock  relations  from 
Equations  1-3  together  with  the  requirement  that  the  entropy  (s)  does  not  decrease.  The 
gas  is  assumed  to  be  adiabatic  and  perfect.  The  total  enthalpy  (ht),  toted  temperature 
(T<),  total  pressure  (pt),  and  total  speed  of  sound  (uj)  are  constant  across  the  shock.  All 
oblique  shock  relations  used  in  this  section  axe  derived  in  Reference  45. 


The  Prandtl-Meyer  function  for  an  isentropic  flow  expansion  is  used  to  determine  the 
flowfield  properties  downstream  of  a  curved  surface  (convex  or  concave)  where  shocks  are 
not  involved.  Locally  two-dimensional  flow  is  assumed.  This  function  cam  be  used  over  a 
shairp  concave  comer  (e.g.  cylinder-flaire  junction),  although  it  is  not  formally  valid  since 
the  flow  is  not  isentropic.  As  a  result,  a  second-order  correction  is  required  and  will  be 
discussed  latter  in  this  section.  For  a  perfect  gas,  the  Prandtl-Meyer  angle,  t/,  through 
which  the  streaun  turns  in  expamding  from  M  =  1  to  a  supersonic  Mach  number,  M,  is 
given  in  Reference  45  ais 


u  =  I  arctan  ^  ~  I)  ~  arctan  y/M'^  —  1 


(16) 


For  first-order  shock-expansion  theory,  we  begin  the  inviscid  flowfield  solution  at  the 
nose  of  the  body.  For  a  nose  half-amgle,  amd  freestream  Mach  number.  Moo,  the  shock 
wave  haJf-angle,  the  Mach  number  behind  the  shock,  M„,  and  the  static  pressure  jump 
across  the  shock,  Pn/Poo,  can  be  determined  from  the  oblique  shock  relations:^* 


cot  6n  =  tan 


h+mi 

2(M^sin*^„  -  1) 


(17) 


2  _  (7  +  1)^M^  sin^  9„  -  4(M2,  sin^  sin^  6n  1) 

"  "  [2^Ml  sin^  -  (7  -  1)][(7  -  IWI  sin»  -h  2] 

Pn  _  27M^  sin^  gn  -  (7  -  I) 

Poo  7  +  1 


(18) 

(19) 


*^E§ftT$,  A.J.,  Stvin,  R.C.,  tnd  Sfneriten,  C.A.,  “The  Genenlieti  Shock- Etpeneion  Method  end  It*  ApfKeotio*  to  Bodie* 
Traoeling  *t  High  Sopereonie  Airipeedi,"  Jotmol  of  Aeronontieol  Science*,  Vol.  tt,  No.  4,  April  19SS,  pp.  tSl-t3S. 

** Spoertoo*,  C.A.,  end  Denni*,  D.H.,  “A  Seeond-Order  Shoek-E*pon*ion  Method  AppKeohle  to  Bodie*  of  Berolutiou  Neor 
Zero  Lift,"  NACA  TN-SSi7,  Jon.  1958. 

*^Ame*  Reoeoreh  Stoff,  “Egnotion*,  Tohle*,  ond  Ckort*  for  Compre**Ule  Flow,"  NACA  Report-I tSS,  NASA  ^me«  Aero- 
noutieol  Lohontorp,  Moffet  Field  CA,  1959. 
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Equation  17  is  solved  iteratively.  The  shock  wave  must  be  attached  to  the  nose  of  the  body 
so  that  0  <  ^  where  S^ax  is  determined  iteratively  from  the  following  equations:^^ 


sin’  6^.^  = 


(7  -t-  l)Ml,  -4  +  yJ(y  +  1)[(7  +  1)M^  +  8(7  -  l)Ml,  +  16] 


cot  SjDMx  =  tan  6, 


( 


47Mi 

(7+l)Ae 


2(M^sin’d, 


rii-) 


Figure  2  shows  the  nose  half-angle  limit,  Smaxt  for  1  ^  -Woo  <  10. 


(20) 

(21) 


A  value  for  »/„  is  obtained  using  Equation  16  and  Mn-  For  all  body  points,  t,  down¬ 
stream  of  the  nose,  Vi  =  6^  —  Si  +  Vn,  where  6i  is  the  local  body  slope.  Since,  Vi  =  f(Mi), 
Ek^uation  16  can  be  used  iteratively  to  obtain  the  value  of  Mi,  the  local  Meudi  nvunber. 
Other  body  surface  flow  properties,  Cp,  (local  pressure  coefficient),  VifVoo  (local  tzmgential 
velocity),  and  Pi/poo  (local  flow  density)  can  be  obtained  from  Equations  13-15. 

The  second-order  theory^^  accounts  for  non-isentropic  and  axisymmetric  effects  using 
the  first-order  results,  denoted  with  the  subscript  s  (e.g.  p,  and  M,). 


Pi  =  Pconci  -  (Pcone,  “  p,,  )exp(-C0) 

(22) 

(  =  {xi  sin  Sn  —  ri  cos  ^„)/(rj  cos  6,-  —  Xi  sin  Si) 

(23) 

,  'TPsMl  sin6A 

2{Ml  -  1  )(pcone.  -  P,  J  V  sin  Sn  ) 

(24) 

A/I  (7  +  1)/2  j 

(28) 

where  Pcon*.  is  the  pressure  on  a  cone  of  half-angle  6i  for  M^o  (see  III.l.c). 


For  flare  afterbodies,  allowance  must  be  made  for  the  non-isentropic  nature  of  the 
flow  thru  the  shock  wave.  Robinson^  derived  a  second-order  correction  for  flzu’es  with  the 
assiunption  that  dpjdx  =  0  at  the  cylinder-flare  junction,  thus  restricting  these  results  to 
bodies  with  long  cylinder  sections. 

The  effect  of  axisymmetric  flow  on  flare  pressure  values,  p,-,  is  computed  using  an 
exponential  curve  fitted  to  the  pressure  gradient  downstream  of  the  cylinder-flare  junction, 
{dpldx)F. 

Pi  =  PcontF  -  (Pconer  "  P»F )  exp(-r?)  (26) 

dp\  /  X  —  xf 

dx)  p  V(Pconer  -  P.f)  COS  6f 
- Z  /  - sin^F  (28) 

\dx)p  \3  —  tanpF/tan(^F  -  ^f)/ 

**/toHnton,  M.L.,  “Beandaty  Lager  Effeeti  in  Sapereonie  Fhw  Over  Cplinder-FUrt  Bediet,”  WRE-Repari-ltSS,  Aaetralian 
Dt/eate  Scientific  Sereicc,  Weapene  Reeearck  Eetailiekmcnt,  Jnip  1974. 
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where 

xp  =  the  location  of  the  cylinder-flare  junction  from  the  nose  of  the  body 
=  pressure  downstream  of  the  flare  shock  (two-dimensional  calcxilation) 

M,f.  =  Mach  number  downstream  of  the  flare  shock  (two-dimensional  calculation) 

Poantf.  =  cone  pressure  using  flare  angle  Sp 
fip  ='arcsin(l/Af«,) 

Bf  =  7P.,M?,/(2(Af?,  -  1)1 

9p  =  shock  wave  half-angle  on  a  wedge  of  half-angle  6p  and  for  Moo 
c.  Cone  Surface  Pressure 

Both  linearized  perturbation  and  second-order  shock-expansion  theories  require  an 
accurate  calculation  procedure  for  the  cone  surface  pressure  at  the  nose  of  the  body. 
Second-order  shock-expansion  theory  requires  cone  pressure  calculations  for  additional 
body  angles.  Due  to  the  frequency  at  which  cone  pressures  must  be  evaluated  for  a  given 
body,  fast  but  accurate  calculation  methods  must  replace  advanced  cone  flow  theories, 
such  as  that  of  Taylor  and  Maccoll.'*^ 

First-order  slender-body  theory  is  accurate  for  cone  half-angles,  6  <  5°  in  supersonic 
flow  {Moo  <  4)  and  prescribes  the  cone  surface  pressure  coefficient  as 


Second-order  slender-body  theory  is  accurate  for  S  <  10®  in  supersonic  flow  {Moo  <  4) 
and  imposes  a  mathematical  limit  on  6  and  Moo  (Note:  coi6/^  >  1  and  cot^^  —  >  0) 

similar  to  linearized  perturbation  theory  (see  Equation  12).  The  cone  surface  pressure  is 
prescribed  as 


where 


+  tan  6  cosh~*(cot  6/^)] 

Linearized  perturbation  theory  could  also  be  used.  However,  slender-body  theory 
(Equ.  29-30)  as  well  as  linearized  perturbation  theory  have  limitations  in  flow  deflection 
angle  and  Mach  number  that  are  imacceptable  for  a  general  projectile  design  code.  Two- 
dimensional  shock  theory  for  surface  pressure  on  a  wedge  is  is  more  general  (see  Figure 
2),  but  does  not  include  the  necessary  axisymmetric  effects  that  are  present  in  cone  flow. 
A  method  that  corrects  the  wedge  theory  for  a  cone  would  be  the  best  choice.  Using 
Equations  17  thru  19  for  the  wedge  theory  and  Reference  47  for  the  cone  theory.  Figure 

Ttflor,  G.L,  Mi  Mmeeoll,  J.W.,  “Tke  Air  Pn«t«rc  »»  •  C«»t  tt  Hi$k  Sfiti,"  Proe.  S«e.  (L»%i»%),  Str. 

A,  VoL  fp.  tli-Slt. 


u  =  — gcosh"^ 

V  =  ayjcxii'^8  — 

a  =  14otan6/[^cot*6  — 
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3  shows  vs.  logio(Moo^)  for  5°  <  6  <  30®  and  2  <  Moo  <  10.  The  data  can 

be  fitted  by  a  fifth-order  polynomial  using  a  least-squares  technique.  Thus  can  be 

determined  from 


d.  Modified-Newtonian  Theory 


The  local  pressure  coefficient  at  a  point  i  on  the  body  can  be  estimated  from  the 
freestream  Mach  number,  Moo,  ratio  of  specific  heats,  7,  and  the  local  body  surface  incli¬ 
nation,  Si,  using^* 

Cp.  =  (31) 


The  body  surface  inclination  angle  is  given  by  sinSi  =  (Ko/  1  loo  1)  •  Other  body  surface 
flow  properties  at  point  i  on  the  body  are  obtained  from  Pn,  amd  Mn  where  subscript 
n  denotes  static  flow  properties  behind  the  projectile  nose  shock. 


Pi  _  C,XPooVI/2)  +  Poo 

poo  \ 


Mi  = 


2 

.7-1 


(  1  +  ^Ml 


(33) 

(34) 


The  Newtonian  model  assumes  that  the  freestream  flow  is  a  streaun  of  particles  in 
rectilinear  motion.  Upon  contact  with  the  body,  this  stream  looses  momentum  normal  to 
the  body  surface  but  retaiins  momentum  tangential  to  the  surface.  These  aissumptions  are 
reasonably  similar  to  the  case  where  the  freestream  Mach  number  is  large,  i.e.  hypersonic 
flows.  Modified-Newtonian  theory  cam  be  applied  to  blunt  bodies.  For  blunt-nose  pro¬ 
jectiles  modified-Newtonian  theory  can  be  combined  with  second-order  shock-expansion 
theory  to  provide  inviscid  flow  surface  profiles  over  the  entire  projectile.  In  contraist,  lin¬ 
earized  perturbation  theory  can  be  applied  solely  to  pointed-nose  projectiles. 


2.  Viscous  Methods 


The  viscous  flowfield  is  considered  to  be  limited  to  projectile  surface  boimdary  layers. 
For  steawiy  two-dimensional/aodsymmetric  flow,  the  governing  botmdary-layer  equations 
describe  the  conservation  of  maiss,  momentum,  and  energy  in  a  compressible  flow. 


djr'^pu)  djr’^pv)  _ 
dx  dy 


_ du  _du 


M  +  L— 

dx  r*  dy 


**  AnitT$tn,  J.D.,  HntT$o%ie  and  High  Tempenttre  G*t  Dfnmmiei,  MeGniu-ffill,  1969. 


(36) 
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_ dH  f-,  ^ 

dz  * i''(‘-p;j“8i;+p;^"'~M. 


_dH  _dH 

pu 


(37) 


where  fc  is  1  for  axisymmetric  flow  and  0  for  twodimensional  flow. 


The  boimdary  conditions  for  Equations  35-37  are  prescribed  as  zero  velocity  on  the 
walls  and  Uniy_*y,  u(i,  y)  =  Ue(i),limy_y.  h(i,y)  =  he(x)  at  the  boimdary  layer  edge. 
The  boimdary  layer  edge  velocity  and  enthalpy,  Ug  and  hg,  are  prescribed  as  well.  The 
reference  velocity  is  the  freestream  value,  Voo-  The  coordinate  x  denotes  the  distance 
along  the  surface  measured  from  the  stagnation  point  or  the  leading  edge  of  the  body.  The 
coordinate  y  is  measured  normal  to  the  body  surface,  which  is  at  a  radius  of  r  from  the 
body  centerline.  The  equations  are  written  for  turbulent  boundary  layer  flow  and  thus  the 
overbar  (u)  represents  a  time- averaged  quantity  zmd  the  turbulent  trzmsport  properties  are 
represented  by  —pu'v'  and  —pv'h'.  Closure  assumptions  must  be  made  for  these  terms. 
The  Boussineq  eddy-diffusivity  concept^*  is  used  for  that  purpose. 

-  pu'v'  =  p€m{du/dy)  (38) 

-  pv'h'  =  p{em/PTtr){dh/dy)  (39) 


The  eddy  viscosity,  is  defined  in  terms  of  an  inner  and  8in  outer  region  within 
the  boundary  layer  and  is  based  on  the  mixing  length,  Lm,  concept.®®  The  concept 
does  not  apply  within  the  sublayer,  therefore  the  calculation  is  started  at  a  sm£dl  distance 
from  the  wall,  yg.  An  intermittency  term,  u)tr,  is  often  introduced  such  that,  uitr  =  0 
for  fully  laminar  flows,  1  for  fully  turbulent  flows,  and  a  function  of  x  for  transitional 
flows.  Transitional  flows  are  not  considered  here.  The  turbulent  Prandtl  number,  Prtr,  is 
a  constant  .9  throughout  the  boundary  layer  which  is  a  good  mean  value  for  air. 

dll 

(€m)inner  =  LI^T  \  -^  \  for  0  <  y  <  yg  (40) 

ay 

(€m)outer  ~  O  |  ^  |  ^tr  for  yc  —  V  —  ye  (^^) 


where, 

a  =  .0168 

S-  =J^{l-ipu}/{p,u,))dy 
Lm  =  xy[l -exp(-y/r)] 

X  =  .4,  the  von  Karmzin  constsmt 
r  =  Van  Driest  damping  length  constant 


For  continuous  boundary-layer  profiles,  (e„)inner  is  applied  until  (Cm)uuieryc  =  (cm)<»uter- 
Cebeci'*®'*®  generalized  the  definition  of  the  Van  Driest  damping  length  constant,  P,  to 
compressible  boundary-layer  flows  with  pressure  gradients,  heat  and  mass  transfer. 


(42) 


**Ce4eei,  T.,  Smith,  A.M.O.,  Analftit  of  T%rh%Unt  Bo%nitrt  Ltteri,  Aetdemie  Prti$,  1974- 

^ Bndekaw,  P.,  Cekect,  T.,  and  Wkiielavj,  J.H.,  Eniineerint  Caletlatitn  Metkodt  lor  Ttritleni  Flovi,  Academic  Prt$t, 
I9SI. 
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where 

N  =  [1  -  11.8{flu,/He){pe/pwyP^Y^^ 

p+  =  (/ieWe)/(PeUr)  dujdx 

Ur  =  {twIPwY^'^ 

F'*’  =  empirical  constant  which  can  be  expressed  as  a  function  of  Reynolds  number. 


Equations  35-37  apply  to  attached  laminar  and  turbulent  botmdary-layer  flows.  Ef¬ 
forts  to  solve  Equations  35-37  with  flow  sepeiration  (i.e.  negative  wall  shear)  are  compli¬ 
cated  both  by  the  need  to  accoimt  for  downstream  conditions  after  separation  and  by  an 
apparent  singularity  at  the  separation  point.  There  is  considerable  analytical  and  numer¬ 
ical  evidence  that  such  a  singularity  exists  when  the  steady  botmdary-layer  equations  are 
solved  in  the  conventional  manner  with  a  specified  boundary  layer  edge  velocity  or  pressure 
distribution.  The  separation  singularity  heis  been  discussed,  for  example  by  Brown  and 
Stewartson.®*  Numerical  integration  of  the  boundary-layer  equations  past  the  point  of 
flow  separation  gives  rise  to  an  instability  in  the  computation.  To  continue  the  solution 
downstream,  it  is  therefore  necessary  to  impose  an  asymptotic  boundary  condition  of  some 
kind  (usually  for  the  pressure)  as  r  — ♦  oo.®^  A  stable,  but  approximate,  downstream  solu¬ 
tion  can  be  obtained  without  such  a  boundary  condition  using  the  FLARE  technique  first 
proposed  by  Flugge-Lotz  And  REyhner;®^  They  neglect  the  '^{du/dx)  term  in  Equation 
36,  which  is  the  cause  of  the  instability,  whenever  u  <  0.  Williams®^  used  the  FLARE 
approximation  in  an  iterative  solution  procedure  called  DUIT.  This  procedure  solved 
the  boundary-layer  equations  using  both  a  standard  and  an  inverse  method  (see  Section 
III.2.b)  and  links  these  solutions  together  by  an  iteration  method.  The  FLARE  approxi¬ 
mation,  used  with  a  standard  solution  of  the  boundary-layer  equations,  accurately  predicts 
the  pressure  at  the  boundary  layer  edge,  but  not  the  skin  friction.  However,  Cebeci®®  and 
Carter®®  have  shown  that  iterative  procedures  employing  FLARE  can  accurately  predict 
the  skin  friction.  FLARE  approximations  have  eJso  been  investigated  by  Kaufman®^  and 
McDonald  amd  Briley.®*  In  the  present  model,  the  FLARE  approximation  is  employed 
in  both  staindard  and  inverse  boundary-layer  solutions  as 


du 

^■^=Cp\u 


du 

dx 


,  C  =  1.0  when  u  >  0  and  C  <  0.2  when  u  <  0 


(43) 


Despite  the  utility  of  the  FLARE  approximation,  the  belief  that  the  boimdary-layer 
equations  should  nevertheless  describe  thin  sepairation  bubbles  persists.®*  This  idea  has 
led  natmrally  to  the  hypothesis  that  the  singularity  results  from  a  failure  to  allow  in  some 

S.N.,  »ud  Siewsrlton,  K,,  “Ltmintr  Siftraiion,"  Annntl  Rtvitw  oj  FUii  Mtektniet,  Vol.  t,  1 969. 

SUwmrtsou,  K.,  ami  WilKtmt,  P.G.,  Vol.  to,  1973,  pf.  96-106. 

*^Fhtit-LoU,  I.,  and  Repkner,  T.A.,  “Tke  Interaction  of  a  Skock  Ware  witk  a  Laminar  Bonndarp  Laptr" 
International  Jonraal  oJ  Non-Linear  Meekaniet,  Vol.  3,  1986,  pp.  173-199. 

^Williamo,  P.G.,  “  A  Reoeroe  Flow  Compntation  In  Tke  Tkeorp  of  Self-Indneed  Separation, "  hectare  Note*  i»  Pkpoiet, 
Vol.  3S,  Proceedinpo  of  tke  4tk  International  Conference  on  Namerieal  Metkodo  i»  Flaid  Dpnamieo,  Sprinier-Verlag,  1979. 

^^Cekeei,  T.,  "Separated  Flow*  and  Tkeir  Repretentaiion  hp  Bonndarp  Laptr  Epnation*,"  Report  ONR-CR  tlS-tSd-t,  Office 
of  Nana!  Reeearck,  Arlington,  VA,  1978. 

** Carter,  J.E.,  "A  New  Bonndarp-Laper  Interaction  Tecknigne  for  Separated  Flow*,"  NASA  TM-76890,  Jane  1976. 

Kaafmann,  K.C.,  an^  Hoffman,  G.H.,  "A  Unified  Approaek  fo  Direct  and  Inveree  Bonndarp  Laper  Solation*," 
International  Joamal  for  Namerieal  Metkod*  in  Flaidt,  Vol.  6,  1966,  pp.  1387-1366. 

^ McDonald,  H.,  and  Brilep,  fV.R.,  “A  Sarvep  of  Recent  Work  on  Interacting  Bonndarp-Laper  Tkeorp  for  Flow  Witk 
Separation,"  Namerieal  and  Pkpeical  Aepect*  of  Aerodpnamic  Flow*  II,  (Cekeei,  ed.)  Springer- Verlag,  1964. 

Brilep,  W.R.,  ani  McDonald,  H.,  "Namerieal  Prediction  of  Ineomprenihle  Separation  Bakkiet, "  Joamal  of  Flaid 
Dpnamiei,  Vol.  89,  Pt.  4,  1979,  pp.  831-696. 


15 


manner  for  interaction  with  the  inviscid  flow.  Since  the  usual  singularity  is  manifested 
by  an  infinite  boundary  layer  displacement  thickness  and  an  infinite  streamwise  gradient 
of  the  wall  shear  at  separation,  efforts  have  been  made  to  compute  inverse  solutions  to 
the  boimdary-layer  equations  by  specifying  a  smooth  displacement  thickness  and  comput¬ 
ing  the  resultant  velocity  distribution.  There  is  extensive  numerical  evidence  that  this 
procedure  removes  the  singularity  at  separation. 

Equations  35-37  can  be  solved  in  a  standard  or  inverse  procedure,  with  each  having 
particular  applications  and  advantages.  In  the  standard  solution,  the  pressure  gradient 
(dp/dx)  imposed  upon  the  boundary  layer  by  the  inviscid  flo^eld  is  specified.  The  numer¬ 
ical  integration  of  the  equations  then  yields  t((x,  y),  r(x,y),  /»(x,y),  p(x,y),  6*,  and  C/ = 
2T^/{pulg^)  where  =  fjidu/dy.  In  the  inverse  solution,  the  boundary  layer  displacement 
thickness,  6*,  or  skin  friction,  C/,  distribiUion  is  specified.  The  numerical  integration  of 
the  equations  then  yields  u(x,  y),  tJ(x,  y),  h(x,  y),  p(x,  y),  and  pe  or  Ug  distributions.  The 
boundary  layer  edge  pressure  and  velocity  are  related  by,  —dp/dx  =  peUedue/dx,  (Euler’s 
equation). 


a.  Standard  Solution  Procedure 

The  boundary-layer  equations  can  be  numerically  integrated  using  the  CS  method 
by  Cebeci  and  Smith. In  the  CS  method,  Equations  35-37  are  transformed  into  nearly 
two-dimensional  form  using  the  Mangier  transformation  and  the  Levy- Lees  transformation 
in  order  to  remove  the  singularity  at  i  =  0  and  stretch  the  coordinates  in  the  x  and  y 
directions.  The  combined  Mangler-Levy-Lees  transformation  is  given  by 

di  -  pePeUt[r{x)l Ltf^dx  and  drf  -  [pUe/i2iY^^][r{x)/Lr]'‘dy  (44) 

where  L,  is  a  reference  length.  The  stream  function,  V’?  is  defined  in  terms  of  p,u,  and  v. 
A  dimensionless  stream  function,  /,  is  related  to  ^  as  follows: 

where  ^  ~ 

Equation  35  is  automatically  satisfied  by  the  definition  of  ip.  Equations  36  and  37  become 
the  transformed  momentum  and  energy  equations  for  two-dimensional  and  axisymmetric 
laminar  and  turbulent  bound2U'y  layers.  For  convenience  the  overbar,  (u),  notation  has 
been  dropped. 

(»/")'+  //"  +  Ale  -(/)>]  =  2f  [{'  ^  («) 

(<..«'  +  a,//")'  +  fg'  =  2i  (47) 

where, 

O' 

/  = 

h  =(l-hr)«A'(l-h4) 
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«m  = 

A  =  (2$/Ue)(ciUe/dO 

C  =  Pe/P 

g  =  hIK 

a,  =  (1  +  r)2*(ii:/Pr)[l  +  6+ (Pr/PrO] 

a,  =(l+r)2''A:(u2/he)[l-(l/Pr)] 

T  =  the  transverse  curvature  parameter,  0  for  zero  yaw 

The  boundary  conditions  and  eddy  viscosity  formulae  are  also  transformed.  Equations 
46  and  47,  for  a  given  pressure  gradient  parameter  A,  contain  only  derivatives  in  the  body- 
normal  direction,  r}.  The  equations  can  therefore  be  solved  using  a  marching  technique 
from  an  initial  solution  at  ^  =  0  to  any  location  downstream  on  the  body.  It  is  exactly 
this  characteristic  that  leads  to  the  instability  for  separated  flows,  discussed  above,  and 
for  which  the  FLARE  approximation  is  applied  to  the  f'df'  fd^  term  in  Equation  46. 
The  initial  solution  (at  ^  =  0)  can  obtained  from  the  compressible  similar  boundary-layer 
equations  for  which  solutions  are  well  known.^® 

The  CS  method  uses  the  Keller  Box  Scheme  to  numerically  integrate  Equations  46  and 
47.  The  scheme  is  accurate  and  efficient  for  solving  laminar  and  turbulent  boundary-layer 
equations.  The  scheme  was  developed  by  Keller,®®  and  the  application  to  boimdary-layer 
equations  was  given  by  Cebeci  and  Smith.^®  One  of  the  basic  features  of  the  scheme  is 
to  write  the  governing  system  of  partial  differential  equations  in  the  form  of  a  first-order 
system.  Thus  derivatives  of  some  quantities  with  respect  normal  variable,  rj,  must  be  in¬ 
troduced  as  new  unknown  functions.  Derivatives  with  respect  to  all  other  (streamwise) 
variables  occur  only  to  first  order,  as  a  consequence  of  the  boundary-layer  approxima¬ 
tions.  With  the  resulting  first-order  system  zind  the  use  of  an  arbitrary  rectangular  grid, 
centered- difference  quotients  and  averages  at  the  midpoints  of  the  grid  cells  yield  second- 
order  accurate  finite-difference  equations.  The  scheme  is  unconditionally  stable,  but  the 
equations  are  implicit  and  highly  nonlinear.  Newton’s  method  is  used  as  a  solution  tech¬ 
nique  along  with  block-tridiagonal  factorization. 


b.  Inverse  Solution  Procedure 


Inverse  solution  of  Equations  46  and  47  has  been  shown  to  be  usefiJ  in  the  numerical 
integration  of  the  bormdary-layer  equations  past  the  point  of  vanishing  wall  shear  stress 
(i.e.  separated  flow)  by  Cebeci.®*  In  the  inverse  solution  procedure,  the  boimdary  layer 
displacement  thickness,  8*,  or  skin  friction,  C/,  distribution  is  specified. 


Keller,  H.B.,  “A  New  Difference  Scheme  Jor  Paniolie  Prohleme,”  Nnmericnl  Sointioni  oj  Pertial  DifferentinI  Etnntieue 
II,  (Bmmhie  ei.)  Aeeiemie  Preis,  1970. 

**  Ceteei,  T.,  “An  Inreree  Botndarf-Lnfer  Method  for  CompreteHle  Lominor  «ni  Tnrhulent  Bonndmrp  Lopere,"  A/AA 
Jotmtl  of  Airenft,  Vol.  IS,  No.  9,  Sept.  1970,  pp.  709-717. 
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The  numerical  integration  of  the  boundzu'y-layer  equations  then  yields  the  boundary  layer 
profiles  as  well  as  the  edge  velocity  (from  which  the  edge  pressure  can  be  computed  using 
Euler’s  equation,  Equation  4). 

To  describe  the  inverse  solution  procedure  of  Cebeci,  assume  (or  solve  for  using  the 
standard  solution  procedure)  the  quantities  /,  /  ,  /  ,  g,  g  ,  A,  and  Ug  at  z  =  z,_i .  At  x  =  x< 
the  solution  of  Equations  46  and  47  for  a  given  S*  is  required.  To  start  the  calculation  it 
is  necessary  to  know  A(xi)  and  Ue(i»)-  ^  value  of  A(xj)  is  chosen  (possibly  extrapolated 
from  upstream  values),  and  a  value  for  Ue(xi)  is  cadculated  from  the  definition  of  A.  Using 
a  centred  finite-difference  for  A(xj_i/2)  terms  of  A(xi)  and  A(xj_i),  Ue(xi)  is  determined 
obtained  using 

Ue(li)  =  -Ue(Xi_i) 

where  Ui  =  Xi-i/2/lxi  -  x._i]  and  A(xi_i/2)  =  [A(a;0  -  A(x,_i)]/2 


A(x._i/2)  +  2a>.- 


(50) 


Once  A(xi)  and  Ue{xi)  are  known,  then  the  stemdard  solution  procedure  can  be  used  to 
solve  for  the  boundary  layer  profiles  and  S’  (denoted  (5j)  at  x  =  x,.  Recall  that  the  desired 
value  for  the  displacement  thickness  is  <5*.  The  function  n  is  then  formed; 


^[A(x.)]  =  s;  -  S’ 


(51) 


A(x,)  is  determined  such  that  7r[A(x,)]  =  0.  This  roctfinding  problem  can  be  efficiently 
solved  using  the  secant  method.®^  For  two  initial  guesses,  A°(ii)  and  A^(ii),  an  iteration 
procedure  can  be  formed  such  that 


”'lA'(xi)|  -  "•[A'*’(''i)] 


(A'(x,)-A-'(x,)) 


(52) 


Equation  52  converges  rapidly  for  most  applications  and  yields  the  values  of  A(x,  )  Eind  thus 
Ue(x,)  for  the  specified  value  of  the  displacement  thickness,  S’{xi).  The  same  procedure 
can  be  used  for  a  specified  C/. 

Cebeci®^  has  shown  that  the  inverse  solution  procedure  yields  good  results  for  nearly- 
separated  (i.e.  negligible  wall  shear)  and  sepeirated  (i.e.  negative  wall  shear)  flows,  when 
compared  to  data  measured  for  decelerating  flows  and  flows  with  strong  adverse  presstire 
gradients.  In  some  cases  the  inverse  solution  method  yielded  results  that  were  superior 
when  compared  to  those  of  the  standard  solution  method. 

Carter®^  has  employed  an  inverse  approach  for  an  incompressible  boundaury-layer  solu¬ 
tion  using  a  different  iteration  procedure  than  described  above.  This  approach  is  adequate 
for  solving  two-dimensional,  as  well  as  three-dimensional,  leiminar  boimdary  layers  that 
were  undergoing  separation  and  reattachment .®® 

*^Audtrien,  D.A.,  Tannekill,  J.C.,  and  PUichtr,  R.H.,  Fl%id  Dtnamiet  and  Hat  Tr€n$fer,  Hcmitpktn 

Puklitkini  Corp.,  McGraw-HiU,  1984. 

**Cekeei,  T.,  Kktlil,  E.E.,  and  WkileUw,  J.H.,  “CalctUtion  of  Separated  Boundtrp-Laptr  FIoiob."  AIAA  Vol.  17, 

No.  It,  Dec.  1979,  pp.  tt9l-tt9t. 

**  Ctritr,  J.E.,  “Invene  Solniione  For  Lamintr  Boondorp-Loper  Flotus  IVitk  Separation  and  heaUaekment,”  NASA  TR 
R-idZ  Noo.  1975. 

Edviarde,  D.E.,  Carter,  J.E.,  and  Smitk,  F.T.,  “AnalpetB  of  Tkree- Dimensional  Separated  Floui  niitk  tke  Bonndarp-Laper 
Eonations.”  AIAA  Jonmat.  Vol  tS,  No.  S,  Marck  1987,  pp.  380-387. 
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3.  Interaction  Method 


For  aerodynamic  design  calculations  it  is  common  to  obtain  the  pressure  distribution 
about  a  body  using  an  inviscid  flowfield  solution.  The  inviscid  solution  provides  the  nec¬ 
essary  boimdary  condition  for  solving  the  boundary-layer  equations.  For  large  freestream 
Reynolds  numbers  and  smooth  body  geometries,  the  presence  of  the  viscous  boundary 
layer  only  slightly  modifies  the  the  inviscid  fiow  over  the  body.  An  improved  inviscid  flow 
solution  can  be  obtained  by  augmenting  the  physical  profile  of  the  body  with  the  boundary 
layer  displacement  thickness.  The  improved  inviscid  edge  velocity  distribution  can  then 
be  used  to  obtain  yet  another  viscous  flow  solution.  In  principle,  this  viscous-inviscid  in¬ 
teraction  procedure  can  be  continued  in  an  iterative  fashion  imtil  the  pressure  distribution 
about  the  body  converges.  In  practice,  however,  severe  underrelaxation  of  the  changes 
from  one  iterative  cycle  to  another  is  often  required  to  achieve  convergence.®^  For  high 
Reynolds  number  flows,  where  the  boundary  layer  is  attached,  the  improvement  in  solu¬ 
tion  accuracy  gained  from  viscous-inviscid  interaction  is  negligibly  small.  In  these  cases 
it  suffices  for  engineering  design  purposes  to  compute  the  inviscid  and  viscous  flows  inde¬ 
pendently,  in  order  to  determine  the  pressure  and  skin  friction  drag,  respectively.  Flows 
which  separate  or  contain  sepaa-ation  bubbles  are  a  noteable  exception.  The  displacement 
effect  caused  by  the  separated  regions  of  the  flow  locally  alter  the  pressure  distribution  in 
a  significant  manner.  A  rapid  thickening  of  the  boundary  layer  under  the  influence  of  an 
adverse  pressure  gradient  (e.g.  near  a  flared  afterbody),  even  without  separation,  can  also 
alter  the  pressure  distribution  to  an  extent  that  a  reasonable  solution  cannot  be  achieved 
without  viscous-inviscid  interaction.  Indeed,  without  accounting  for  displacement  effects, 
the  solution  obtained  may  predict  flow  separation  when  the  flow  does  not  separate  at  all. 

The  essential  elements  of  a  viscous-inviscid  interaction  calculation  procedure  are:  1)  a 
method  for  obtaiining  the  inviscid  flowfield  about  a  body  of  arbitrary  cross  section  in  order 
to  swrcount  for  viscous  flow  displacement  effects,  2)  a  technique  for  obtaining  a  solution 
to  the  boundary-layer  equations  suitable  for  the  problem  at  hand  -  an  inverse  boimdary- 
layer  procedure  would  be  appropriate  for  a  flow  which  may  separate,  and  3)  a  procedure 
for  relating  the  inviscid  and  viscous  flow  solutions  in  a  manner  which  will  drive  the  changes 
from  one  iteration  cycle  to  the  next  toward  zero.  The  techniques  required  for  items  1  and 
2  have  been  described  in  the  previous  sections.  The  interaction  procedure  is  described  in 
this  section. 

A  number  of  viscous-inviscid  interaction  schemes  have  been  proposed.  Briley  and 
McDonald*®  used  a  time-dependent  interaction  scheme  for  incompressible  flows  with  thin 
separation  bubbles.  They  also  surveyed  a  wide  range  of  interaction  schemes.  Carter*® 
developed  an  interaction  scheme  for  transonic  and  supersonic  separated  flows  that  was 
xised  for  incompressible  separated  flows  by  Kwon  and  Fletcher.^  This  later  scheme  is 
employed  in  the  present  model. 

The  interaction  calculation  proceeds  in  the  following  way.  First,  the  inviscid  edge 
velocity  and  the  boundary  layer  flow  are  computed  up  to  the  beginning  of  the  specified 
interaction  region  {x  =  xi).  A  standard  solution  to  the  boimdary-layer  equations  is  em¬ 
ployed.  These  two  solutions  do  not  change.  Next,  an  initial  botmdary  layer  displacement 
thickness,  ^*(x),  distribution  is  chosen  over  the  interaction  region,  ii  <  x  <  Xj.  This 
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mitial  guess  is  purely  arbitrary,  but  it  should  match  the  ^*(x)  computed  by  the  standard 
method  at  x  =  Xi .  The  boimdary-layer  solution  is  obtained  by  an  inverse  method  using 
this  S*(x)  as  a  boundary  condition,  yielding  the  edge  velocity  distribution,  Ue^(x).  Now 
the  inviscid  solution  is  used  to  obtain  Ue,tn«(x)  using  the  body  profile  augmented  by  the 
displacement  thickness,  ^*(x).  The  Ue(x)  from  the  two  calculations,  boundary-layer  and 
invisdd)  will  not  agree  \mtil  convergence  has  been  achieved. 

The  difference  between  Ue(x)  calctdated  for  inviscid  and  boimdary-layer  flows  can  be 
used  as  a  potential  to  calculate  an  improved  distribution  for  6*{x).  A  suitable  scheme 
follows  from  the  observation  that  the  volume  flow  rate  per  unit  width  in  the  boundary 
layer  tends  to  be  preserved,  i.e.  ex  constant.  This  implies  that  a  local  decrease  in 
Ue(x)  (associated  with  a  more  adverse  pressure  gradient)  causes  an  increase  in  ^*(x)  and 
a  local  increase  in  Ue(x)  (associated  with  a  more  favorable  pressiure  gradient)  causes  a 
decrease  in  ^*(x).  This  concept  is  used  to  compute  the  new  distribution  of  S*(x)  using 

(4^)  (53) 

\  ^e,inv  ) 

where  i  denotes  the  iteration  level.  At  convergence  Ue,6/  =  Ue.mv;  thus  Equation  53  repre¬ 
sents  an  identity  thereby  having  no  effect  on  the  final  solution.  Equation  53  only  serves 
as  a  basis  for  correcting  8'  between  iterative  passes.  No  formal  justification  for  the  use 
of  Equation  53  is  required  so  long  as  the  iterative  process  converges.  However,  Carter*^ 
has  derived  a  formal  justification  of  Equation  53  based  on  the  von  Karman  momentum 
integral  for  boundary-layer  flows.  Kwon  and  Pletcher**  have  found  that  convergence  can 
be  accelerated  by  augmenting  Equation  53  with  successive  over-relaxation  (SOR): 


(6‘)‘+‘ 


=  RF 


+  {l-RF){6y 


(54) 


where  RF  is  the  relaxation  parameter.  Kwon  and  Fletcher^*  fotmd  that  1.2  <  RF  <1.8 
produced  adequate  convergence  acceleration. 


The  viscous-inviscid  interaction  procedure  is  completed  by  making  successive  passes 
first  through  the  inverse  boimdary-layer  scheme,  then  through  the  inviscid  flow  scheme 
with  6*  being  computed  by  Equation  53  and  54  prior  to  each  boundary-layer  calculation. 
Convergence  is  said  to  have  occurred  when 


— ”*’*”*'  '  <  the  prescribed  tolerance 

We,»nv 


(56) 


4.  Small  Yaw  Procedure 

The  calculation  of  the  surface  pressures  on  a  body  of  revolution  is  considerably  more 
complicated  when  the  body  is  inclined  to  the  flow.  When  the  body  is  not  inclined,  the 
meridian  lines  on  the  surface  are  coincident  with  the  streamlines,  and  the  analysis  used 
for  the  nonlifting  body  is  adequate.  At  incidence,  the  streamlines  deviate  considerably 
from  the  meridian  lines  because  of  the  crossflow  component  of  the  flow  on  the  surface.  If 
the  shape  of  the  streamlines  can  be  approximated  and  this  shape  used  to  generate  a  body 
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of  revolution,  it  is  reasonable  to  assume  that  the  analysis  for  the  nonlifting  body  can  be 
used  to  provide  the  pressure  variation  along  a  meridian  of  this  body  of  revolution  and, 
therefore,  along  the  streamline  itself. 

The  approach  taken  in  this  report,  after  Jackson  et.al.,^  is  to  obtain  the  coordinates 
of  the  streamlines  in  the  vertical  plane  {4  =  dO*’)  for  the  body  of  revolution  at  angle  of 
attack,  a,  by  transforming  the  body-axis  coordinates  (x,  y)  into  the  wind-axis  coordinates 
(^^wiadii/wiiHi)*  the  wind-axis  system,  the  circumferential  angle,  is  transformed  to  il 
where  ^  =  0®  and  (1  =  —90®  correspond  to  the  windward  side  and  where  4  =  180®  and 
ft  =  -1-90®  correspond  to  the  leeward  side;  thus  ft  =  ^  —  90®.  This  transformation  results 
in  a  rotation  of  the  coordinate  system.  If  the  body  has  a  blunt  nose  cap  of  radius  r„, 
the  center  of  rotation  is  taken  to  be  the  center  of  the  cap.  The  streamline  for  ft  =  0®  is 
approximated  by  the  body-axis  coordinates  of  a  meridian  for  the  true  body  shape.  These 
asstunptions  imply  that  the  coordinates  of  the  streamlines  on  the  body  at  angle  of  attack 
vary  from  body-axis  coordinates  at  ft  =  0®  to  wind-axis  coordinates  at  ft  =  ±90®.  This 
variation  is  represented  for  —90®  <  ft  <  90°  by  the  expressions: 

/r„(cosa  -  1)  +  X  +  (ycoso  -  (x  -  r„)sinQ;)sina  \  >  /ccx 

a^wind  - - X  smft  +  X  (56) 

\  cosa  J 

ywind  =  (ycoso  -  (x  -  r„)sinQ  -  y)sinft  +  y  (57) 

Using  Equations  56  and  57,  the  coordinates  of  the  streamlines  in  the  wind-axis  system 
are  obtained,  and  these  coordinates  are  used  to  generate  bodies  of  revolution  for  each  cir¬ 
cumferential  angle,  ft,  of  interest.  Estimates  of  pressure  distribution  along  the  body-axis 
meridians  are  obtained  by  applying  the  ansJysis  used  for  the  nonlifting  body  to  equivalent 
bodies  generated  from  the  streamlines.  The  use  of  this  concept  involves  two  basic  limita¬ 
tions:  (1)  the  angle  of  attack,  a,  is  small  enough  for  the  stagnation  point  to  remain  on 
the  nose  cap  (if  the  body  is  blunt)  and  (2)  the  angle  of  attack  is  small  enough  for  the 
coordinates  of  the  equivalent  body  for  ft  =  90°  (leeward  side)  to  remain  positive. 

Such  a  system  of  transformations  result  in  engineering  estimates  of  the  streamline 
shape  as  well  as  the  pressure  variation  along  the  meridians.  An  examination  of  these 
transformations  (Equations  56  and  57)  shows  the  incidence  approaches  zero  (a  — »  0®)  as 
the  approximate  streamline  shapes  approach  the  body  meridian  shapes  (i.e.  true  streamline 
shapes  at  ft  =  0®,  ^  =  90®).  This  approach  should,  therefore,  be  valid  for  the  case  where 
the  body  incidence  is  small  and  crossflow  on  the  projectile  is  negligible. 


IV.  HEAT  TRANSFER  ANALYSIS 


Hypersonic  projectiles  which  attain  speeds  of  2  to  3  km/s  at  sea  level  (Mach  number 
6  to  9)  are  exposed  to  extreme  thermal  and  aerodynamic  loads.  The  stagnation  point  heat 
flux  can  exceed  10  kW/cm*  while  stagnation  pressures  range  from  40  to  100  atms.  These 

**Jtck»on,  C.M.,  S«ivfer,  W.C.,  tni  Smith,  R.S.,  “A  Mttkoi  For  Determinin/  Str/aee  Prtioorti  om  Bl%%t  Boiito  oj 
Beoohtion  *i  Smoll  An0lei  oJ  AtUei  in  Snpertonie  Flow,  ”  NASA  TN  D-4SBS,  Nooemher  IStS. 
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instantaneous  loads  can  be  considerably  larger  than  those  encountered  by  reentry  vehicles. 
However,  the  residence  time  in  this  environment  is  much  less  than  the  corresponding 
time  of  maximum  aerobraking  of  reentry  vehicles.  Consequently,  an  accurate  prediction 
of  the  thermal  loads  is  crucial  in  determining  an  optimal  configuration  and  material  for 
hypervelocity  projectiles.  A  poorly  conceived  design  could  result  in  substantial  material 
loss  and  projectile  deceleration.  This  subsequent  degradation  in  kinetic  energy  could  defeat 
the  ptirpose  of  the  design.  Overall,  the  design  objective  is  to  simultaneously  reduce  the 
aerod}rnamic  drag  and  the  surface  heat  fitix. 

In  general,  the  aerothermodynamic  environment  is  strongly  dependent  upon  the  pro¬ 
jectile  shape,  Mach  number,  and  the  Reynolds  number  (and  Knudsen  number  for  high- 
altitude  fiight).  The  surface  temperature  of  the  projectile,  which  also  effects  the  fiowfield, 
is  a  function  of  the  siuface  heat  fiux,  geometry,  and  physical  properties.  Certain  struc¬ 
tural  effects  result  from  this  rapid  transient  heating.  In  a  precise  analysis,  the  heating 
and  structural  problems  are  coupled,  and  it  is  only  with  certain  approximations  that  the 
two  can  be  treated  separately.  In  this  report  structural  heat  conduction  is  not  addressed. 
The  aerodynamic  heat  transfer  to  the  projectile  surface  is  calculated  assuming  a  constant 
projectile  wall  surface  temperature. 

Computational  procedures  for  aerothermodynamic  prediction  of  high-speed  projectiles 
may  be  categorized  as:  reliable  semi-empirical  methods  based  on  benchmark  solutions 
and  experimental  data  or  advanced  numerical  schemes  (e.g.  real-g£is  CFD  codes)  for  the 
solution  of  the  entire  fiowfield.  The  approximate  methods  yield  results  with  only  a  fraction 
of  the  computational  time  required  for  the  numerical  schemes.  In  addition,  the  code 
development  based  on  these  approximate  methods  is  straightforward,  while  the  numerical 
schemes  usually  require  substantial  code  development,  input  setup,  and  grid  generation. 
Approximate  algebraic  equations  have  been  developed  which  accurately  predict  stagnation 
point  and  laminar  or  turbulent  heat  transfer  distributions  for  axisymmetric  bodies.  These 
equations  assume  that  the  fiowfield  is  a  continuum  and  chemical  nonequilibrium  effects 
are  insignificant.  Projectiles  usually  achieve  significant  heating  rates  at  sufiiciently  low 
altitudes  where  these  zissumptions  axe  reasonable.  When  the  geometry  and  the  subsequent 
fiowfield  are  more  complicated,  the  numerical  schemes  will  be  required. 

Before  presenting  the  formulae  for  heat  transfer  prediction  used  in  this  report  it  is 
appropriate  to  review  some  basic  tenets  of  aerodynamic  beat  transfer.  The  phenomenon 
of  heat  transfer  between  a  fiowing  fluid  and  a  body  immersed  in  the  fiuid  is  a  classical 
engineering  problem.  For  the  purposes  of  this  analysis  aerodynamic  heating  is  defined 
as  the  heating  of  a  body  due  to  the  high  speed  flow  of  a  gas  over  the  body,  including 
negative  heating  (cooling),  of  the  body.  The  basic  principles  of  aerodynamic  heating  are 
contained  in  boundary  layer  theory.  In  the  boimdary  layer,  kinetic  energy  is  transformed 
into  thermal  energy  by  frictional  forces  and  compression  forces.  The  temperature  of  the  gas 
in  the  thin  boundary  layer  is  raised,  suid  the  amoimt  of  heat  transferred  is  proportional  to 
the  difference  between  the  body  surface  temperatiure  and  some  representative  temperature, 
Tr,  of  the  hot  gas.  Denoting  gw  sa  the  heat  transfer  to  the  wall  per  unit  area  and  per  unit 
time  from  a  fluid  of  temperature  Tr  to  a  body  of  surface  temperatmre  Tw,  then  Newton’s 
law  relates  gw  to  a  heat  transfer  coeflScient,  Ch,  as 


=  CH(Tr  -  Tw) 
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The  gas  temperature  Tr  is  called  the  recovery  temperattire  and  is  defined  from  the  fireestream 
temperature,  Mach  number,  Prandtl  number,  and  ratio  of  specific  heats  using^^ 

r,  =  T„  +  Pr(g)=T„(l  +  Pr(l^)wi)  (59) 

At  the  stagnation  point  on  the  body,  Pr  =  1,  and  the  recovery  temperature  is  equivalent  to 
the  total  temperature,  Tt,  of  the  fluid.  The  determination  of  Ck  and  thus  is  complicated 
due  to  the  dependence  on  T^.  In  general,  the  wall  temperature  is  based  on  solving  heat 
conduction  equations  within  the  body  based  on  values  of  as  input.  However,  q^,  depends 
on  Tu,.  Thus,  the  heat  conduction  equation,  which  is  linear,  has  a  nonlinear  boundary 
condition.  However,  the  dependence  of  on  T^/Tr  is  not  very  strong^^.  Thus,  for  a 
first  approximation,  q^  may  be  considered  to  be  independent  of  (or  Ty,  =  constant) 
for  short  time  of  flight  simulations.  More  accurate  results  could  be  achieved  using  a  heat 
transfer /conduction  interaction  method. 

In  this  report,  separate  algebraic  formulae  are  used  to  compute  the  heat  transfer  rates 
at  the  projectile  stagnation  point  and  along  the  projectile  surface  for  laminar  and  turbulent 
boundeuy  layers.  The  boundary  layer  properties  required  for  these  formulae  are  computed 
using  the  INTERACT  code  described  in  section  III. 

Heat  transfer  rates  at  the  stagnation  point  on  a  blunt-nose  projectile  are  calculated 
using  correlation  formulae  for  a  gas  in  chemical  equilibrium.  These  formulae  are  given  by 
Fancett^  based  on  the  work  of  Fay  and  Kemp.^  In  this  formula  the  stagnation  point 
heat  transfer  rate,  in  kW/m^,  is  computed  from  the  freestream  velocity,  Voo  in  m/s,  the 
norm^d  shock  stagnation  pressure,  (p«)n  in  atm,  and  the  blunt-nose  radius,  r„  in  cm.  The 
expression  is  formally  valid  for  a  constant  projectile  wall  temperature  of  300K. 

=  2.428x10-'  (VU)"'”*'  ((p,)„/r„)''=  (60) 

where,  (p,)„  = 

Rubesin  and  Johnson^^  and  Young  emd  Jzmssen^  have  shown  that  if  the  values  of  the 
density  and  viscosity  in  a  boundary  layer  are  calculated  with  relationships  which  hold  for 
constant  property  fluids  at  2m  intermediate  “reference  enthalpy” ,  defined  by  an  empirical 
expression,  then  they  accurately  describe  the  actual  conditions  in  a  supersonic  boundary 
layer.  Eckert^®  used  this  “reference  enthalpy”  technique  to  calculate  heat  transfer  rates  in 
a  zero  pressure  gradient  flow.  He  used  the  empirical  expression  for  “reference  enthalpy” 
temperature,  T*,  suggested  by  Young  and  Janssen^  which  involves  the  body  wall  (T^), 
freestream  (Too),  and  the  recovery  (T,)  temperatures  in  Kelvin. 

r  =  Too  -h  — +  .22{Tr  -  Too)  (61) 

Fancett®*  derived  an  expression  for  the  laminar  and  turbulent  heating  rate  distributions 
over  a  wall,  in  W/m*,  using  these  techniques.  For  a  laminar  boimdary  layer: 

g,„  =  3.32  yT  (Pr*)-"/^  Cp  {Tr-T,,)  (62) 

o 
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For  a  tiirbvilent  boundary  layer: 


=  1.76  (log,„Re;)-"“  c,  (T,  -  T„)  (63) 

where  a  is  the  body  surface  distance  from  the  stagnation  point  in  m,  Cp  is  the  specific  heat 
at  constant  pressure  in  J /(kg-K)  for  air  at  fr^stream  conditions,  /i*  is  the  local  viscosity 
in  kg/(m-s)  based  on  Sutherland’s  law  and  T*,  Re*  is  the  local  Reynolds  number  based  on 
the  local  density,  viscosity  (using  T"),  and  <r,  and  Pr*  is  the  local  Prandtl  number  based 
on  /i*  and  Cp. 

The  specific  heat  at  constant  pressure  used  in  Equations  62  and  63  can  be  formulated 
to  include  real  gas  effects  in  the  form  of  a  temperature  dependence.'*^  For  a  diatomic  gas, 
such  as  O2: 

where,  =  2230K  for  O2,  T  is  the  gas  temperature  in  K,  and  is  the  specific  gas  constant 
in  J/kg/K.  Note,  for  T  -C  ^„,Cp  — ►  (7/2)3i,  and  T  >•  ^„,Cp  — »  (9/2)3i. 


The  Stanton  number  distribution  over  the  projectile  surface  is  computed  from  the  heat 
transfer  rate  distribution.  The  Stanton  number  based  on  boundary  layer  edge  conditions 
and  Stanton  number  bzised  on  freestream  conditions  are  defined  as 


St 


Qw 


ptUgCpiTr  -  Tu,) 

Pe«e 


Stoo  =  Ste 


PooUc 


(65) 

(66) 


V.  CODE  VALIDATION 


1.  Flow  Profiles 

In  this  section  pressure  coefficient  distributions  over  projectile  surfaces  as  predicted 
using  the  INTERACT  code  are  compared  to  measurements  on  wind  tiinnel  models.  These 
models  were  instrumented  with  surface  pressure  taps  and  transducers  with  a  reported 
accuracy  of  about  2%  of  the  freestream  pressure.®^**®  Since  the  INTERACT  code  has  been 
setup  for  solely  axisymmetric  geometries,  fin-stabilized  projectiles  are  not  addressed  in  this 
report.  Instead,  long-rod  projectiles  that  are  flare-stabilized  are  considered.  Preliminary 
design  investigations  for  hypersonic  projectiles  have  also  considered  flare-stabilization  over 
the  use  of  fins,  primarily  due  to  the  enhanced  resistance  to  thermal  dsunage  in  flight. 

Figure  4  shows  the  pressure  coefficient  distribution  along  a  5-degree  cone,  cylinder 
midsection,  and  4- degree  flare  afterbody  tested  in  a  wind  tunnel  at  Mach  =  2  and  a 
Reynolds  nxunber,  based  on  total  body  length,  of  7.6  million  (Re  based  on  body  diameter 

WttKin§ton,  W.D.,  tni  Hnmpkrtt,  J.A.,  “Prt$nrt  Mem$%Ttment§  on  Fo%t  Cone-Ctlinder-FItn  CeHfigurtiient  at  Safer- 
eenie  Sfeede,”  U  S.  Army  Mieiile  Command  Rtfort  No.  RD-TM-S9-11,  tO  Oct.  1969. 

**GTaf,  J.D.,  “Laminar  Boandarf-Laftr  Scfaration  on  Flared  Bodice  at  Snpereonie  onrf  Hffereonie  Sfeed$,“  U.S.  Air  Force 
Sfoteme  Command  AEDC-TR-64-t77,  Jan.  1965. 
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of  .5  million).^^  For  free  flight  at  sea  level  the  Re  would  be  15  million  for  Mach  =  2.  The 
model  was  at  zero  yaw.  The  cone  section  extends  to  x/D  =  5.6,  the  cylinder  section  extends 
to  xfD  =  14.2,  and  the  flare  afterbody  extends  to  the  total  body  length  of  about  16.2 
calibers.  The  ntunerical  prediction  using  the  INTERACT  code  is  shown  with  an  inviscid 
(Re  =  oo)  calculation.  The  prediction  that  couples  the  viscous  effects  with  the  inviscid 
flow  yields  the  best  comparison  with  the  measured  data.  Note  that  both  calculations 
reproduce  the  second-order  effects  over  the  cylinder  and  flare  sections.  Note  also  that 
the  fully-coupled  solution  indicates  a  small  degree  of  boimdary-layer-induced  upstream 
influence  at  the  cylinder-flare  junction. 

Figure  5  shows  the  pressure  coefficient  distributions  along  the  same  configuration  as 
Figure  4  but  at  Mach  munbers  of  2  zmd  4  (Reynolds  number  held  constant).  Both  the 
measmed  data®^  and  the  INTERACT  calculations  indicate  that  the  stabilizing  effect  of 
the  flared  afterbody  (caused  by  pressure  levels  that  are  higher  than  the  cone  pressure)  is 
reduced  with  increased  Mach  number.  Agreement  between  the  code  and  data  is  slightly 
better  for  Mach  =  2  at  a  constant  Re.  For  these  wind  tunnel  tests®^  the  Mach  number 
is  increased  by  increasing  the  supply  (totzd)  pressure  (i.e.  static  pressure  reduced),  while 
holding  the  total  temperature  constant,  thus  resulting  in  a  constant  Re  value.  Whereas 
for  free  flight,  an  increase  in  Mach  number  from  2  to  4  results  in  Ein  increase  in  Re  from  15 
to  30  million.  This  effect  on  the  boundary-layer  thickness  in  an  adverse  pressure  gradient 
can  be  overpredicted  by  the  turbulence  model,  especially  on  the  flare. 

Figures  6  and  7  show  the  pressure  coefficient  distributions  along  projectiles  with  6- 
degree  and  8-degree  flare  afterbodies  for  Mach  numbers  2  thru  5  and  Re  =  7.6  million.®^ 
The  numerical  predictions  using  the  INTERACT  code  reproduce  pressure  levels  Euid  trends 
well.  For  a  flare  angle  of  8-degrees  and  Mach  number  5  (Figure  7),  the  boimdary  layer  over 
the  flare  thickens  rapidly  smd  reduces  flsure  pressures.  The  calculations  do  not  adequately 
simulate  this  effect.  Investigation  reveals  that  steep  flare  angles  for  low  Reynolds  number 
(less  than  10  million)  flows  repeatedly  demonstrate  this  effect  in  wind  tunnel  experiments. 
Although  this  effect  is  included  in  the  code  through  the  viscous-inviscid  coupling  (see 
discussion  below),  the  magnitude  of  the  predicted  effect  is  smaller  them  measured.  However, 
the  flare  pressure  level  at  the  end  of  the  body  is  reproduced  by  the  code. 

This  phenomena  is  further  demonstrated  in  Figure  8.  This  figure  shows  the  pressure 
distribution  measured  in  a  wind  tunnel”  along  a  20-degree  flare  with  Re  =  18  million. 
Two  sets  of  data  are  shown.  For  Moo  =  3,  the  cylinder  section  of  the  body  was  shorter 
and  thus  the  flare  begins  at  xfD  =  5.8,  while  the  Moo  =  8.1  flare  begins  at  x/D  =  6.8. 
The  flares  for  both  models  are  of  equal  length.  This  geometry  represents  a  severe  test  of 
the  INTERACT  code,  since  a  flare  angle  of  this  magnitude  causes  rapid  boimdary  layer 
thickening  immediately  after  the  cylinder-flare  junction.  This  effect  may  not  be  reproduca- 
ble  using  the  current  turbulence  model.  However,  as  an  illustration  of  the  viscous-inviscid 
coupling,  these  cases  were  computed.  For  the  Moo  =  3  case  an  Euler/boundary-layer  so¬ 
lution  and  the  coupled  viscous-inviscid  solution  are  compared  to  the  measured  data.  The 
trend  of  reduced  pressure  over  the  inviscid  result  downstream  of  the  cylinder-flare  junction 
is  modeled  by  the  coupled  solution.  The  overall  computed  pressure  level  on  the  flare  is 
10%  smaller  than  measured.  This  is  produced  by  the  cylinder  flow  that  is  not  fully  recom¬ 
pressed  (i.e.  Cp  <  0)  at  the  flare  junction,  as  required  by  the  second-order  shock-expansion 
inviscid  theory  (section  Ill.l.b).  For  the  Moo  =  8.1  case,  the  coupled  INTERACT  solution 
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compares  well  with  the  measured  data.  Again,  the  magnitude  of  the  initial  pressure  level 
reduction  is  overpredicted,  but  the  overall  flare  pressure  is  reproduced.  The  cylinder  flow 
for  this  cjise  was  nearly  recompressed  (i.e.  Cp  ~  0)  as  required  by  the  theory. 

Figure  9  shows  the  pressure  coefficient  distribution  measured  for  a  6*degree  flared  body 
at  Mach  =  2,  Re  =  7.6  million,  and  one-degree  of  yaw.®^  Computations  using  INTERACT 
for  the  windward  and  leeward  sides  =  0®,^  =  180®)  are  com;  ?*red  to  the  data.  Cone 
and  cylinder  section  presstire  levels  for  the  windward  side  are  slightly  overpredicted  and 
underpredicted  for  the  leeward  side.  Pressure  levels  on  the  flare  afterbody  are  reproduced 
by  the  code.  Recall  that  the  small  yaw  procedure  used  in  INTERACT  (section  III.4) 
assumes  no  cross  flow  over  the  body.  Small  cross  flow  effects  on  the  surface  pressure  are 
illustrated  on  the  cone  and  cylinder  sections.  The  effect  of  these  pressure  differences  on 
the  aerodynsimic  stability  predictions  is  addressed  in  the  following  subsection. 

Overall,  the  ability  of  the  INTERACT  code  to  reproduce  body  surface  pressure  coef¬ 
ficient  distributions  heis  been  demonstrated.  The  wind  tunnel  measurements  available  for 
comparison  were  for  low  Re  values  of  7  to  18  milUon  (relative  to  corresponding  free  flight 
Re  values  of  15  to  60  million)  and  in  some  cases  steep  flare  angles  (>  6®).  The  trend  of 
reduced  pressure  (relative  to  the  inviscid  vzdue)  downstream  of  the  cylinder-flare  junction 
is  simulated  by  the  viscous-inviscid  interaction  scheme.  With  increasing  Mach  number, 
reduced  flare  pressure  levels  are  observed  and  computed.  It  is  importamt  to  remember  that 
these  numerical  simulations  were  performed  in  under  one  minute  of  CPU  time  on  a  VAX 
8600  minicomputer. 


2.  Aerodynamic  Stability 


The  design  of  projectile  configurations  with  sufficient  aerodynamic  stability  to  sustadn 
small-yaw  flight  is  one  of  the  most  essential  requirements  of  the  initial  design  cycle.  Non¬ 
spinning  (or  slowly  rolling)  projectiles  must  be  designed  to  maintain  static  stability  over 
the  entire  Mach  number  range  in  flight.  Therefore,  the  engineering  design  code  must 
reliably  predict  the  location  of  the  center  of  pressure,  Xcp,  as  a  function  of  Mach  number. 
The  INTERACT  code  predicts  surface  pressure  and  skin  friction  distributions  over  the 
projectile  as  a  function  of  the  circumferenti2il  body  angle,  <f>,  for  a  given  Moo  ^d  a.  The 
coefficients  of  pressure  drag  (C/j^),  normal  force  (C/v),  and  pitching  moment  (Cm)  are  then 
computed  from  integrations  of  surface  pressure  using 


Cd,  = 
Cs  = 


Axial  Force  1 

Normal  Force  1 


Cm  = 


qoonD‘^/A 
Pitching  Moment 


/L  rtn 

J  CpT  tan  66L<i>dx 


-1 


f  f  Cprx  cos  ^d^dx 
Jo  Jo 


(67) 

(68) 

(69) 


QooDttD'^ 

Skin  friction  drag  coefficient  can  also  be  computed.  The  moment  coefficient  is  referenced  to 
the  nose  of  the  projectile  auid  is  defined  using  pitch  up  as  positive.  The  center  of  pressure. 
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(70) 


also  referenced  to  the  nose  of  the  projectile,  is  computed  from  Cn  and  Cm  using 

_-Cm 
Cn 

Center  of  pressure  data  has  been  measured  in  a  wind  tunnel  for  a  9.5‘’-cone,  cylinder 
midsection,  and  5‘’-flare  configuration  at  Mach  numbers  3  to  6.3  and  small  yaw.^  The 
Reynolds  number  (based  on  total  length)  varied  as  a  function  of  the  Mach  number  between 
1.2  and  6.8  million.  These  values  of  wind  tunnel  Re  give  rise  to  unusually  thick  boxmdary 
layers  as  observed  in  the  test.^^  Table  2  shows  the  measured  and  computed  Xcp  values 
(referenced  to  the  nose  of  the  projectile)  along  with  the  relative  errors  in  the  predictions. 
The  INTERACT  code  and  the  NSWC-AP  code^  are  used  to  supply  predictions.  The 
codes  were  run  with  a  =  1.0°  whereas  the  measurements  were  taken  at  an  tinreported 
angle  of  attack.  However,  center  of  pressure  is  usually  insensitive  to  the  value  of  a  for 
small  yaw. 


T^ble  2.  VARIATION  OF  CENTER  OF  PRESSURE  WITH  MACH  NUMBER 


Comparison  made  for  9.5"  Cone  -  Cylinder  -  5 

Xcp  (cal)  from  the  nose 

"  Flare  and  o  =  1.0® 

Moo 

Re(L) 

Measured 

NSWC  Code 

%  error 

INTERACT  Code 

%  error 

3.0 

6.8x10^ 

4.71 

4.37 

-7.2 

5.21 

+7.2 

4.2 

5.8x10® 

5.05 

4.55 

-9.9 

5.25 

-1-3.9 

5.0 

3.1x10® 

5.77 

4.72 

-18.2 

5.30 

-8.1 

6.3 

1.2x10® 

5.39 

4.87 

-9.6 

5.08 

-5.8 

ave  ; 

-11.2 

±6.25 

The  NSWC-AP  code  averages  11.2%  error  in  Xcp  for  this  configuration,  underpre¬ 
dicting  the  measured  values.  The  INTERACT  code  averages  ±6.25%  error  in  Xcp  values. 
The  trend  in  measured  center  of  pressme  indicates  that  Xcp  moves  rearward  as  the  Mach 
number  increases  with  an  shift  forward  on  the  projectile  for  Moo  =  6.3.  The  forward  shift 
is  due  to  the  decrease  in  Re  and  thickening  of  the  boimdary  layer  at  the  cyUnder-fiare 
junction  thus  decreasing  flare  effectiveness.  The  INTERACT  code  predicts  the  onset  of 
this  phenomena.  The  INTERACT  code  was  run  using  linearized  perturbation  inviscid 
theory  with  a  coupled  finite-difference  boundary  layer  solution.  The  NSWC-AP  code  uses 
second-order  shock-expansion  inviscid  theory  with  a  de-coupled  empirical  boundary  layer 
solution.  As  noted  in  section  Ill.l.a,  linearized  perturbation  theory  more  closely  approxi¬ 
mates  the  Euler  equations.  Investigations  using  the  INTERACT  code  reveal  that  the  loss 
of  flare  effectiveness  at  Mach  =  6.3  is  a  direct  result  of  viscous-inviscid  interaction  at  the 
cylinder-flare  jimction. 

This  analysis  of  predicted  aerodynamic  stabihty  represents  one  particular  projectile 
configuration.  The  INTERACT  code  should  be  checked  with  other  wind  tunnel  and  free 

**Deunu,  D.H.,  "Tkt  af  Bonnitrf-Ltftr  Stptrtiion  Over  Boiiti  of  RenoUtion  Witk  Conietl  Tail  FUrtt,”  NACA 
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flight  measurements  of  center  of  pressure  as  these  data  become  available.  The  favorable 
comparison  with  measured  data  for  a  5°  flared  projectile  is  encouraging  especially  since 
calctilations  were  performed  interactively  on  a  VAX  8600  minicomputer. 


3.  Heat  Transfer  Profiles 

In  this  section  wall  heat  transfer  rate,  and  Stanton  number,  St,  predicted  us¬ 
ing  the  INTERACT  code  are  compared  to  measurements  on  wina  tunnel  models.  Three 
configurations  are  examined:  a  flat  plate  with  sharp  leading-edge,  a  blimt  cone,  and  a  cone- 
cylinder-flare.  The  surface  heat  transfer  data  referenced  in  this  report,  were  measured  on 
wind  tunnel  models  using  either  temperature  sensitive  paint  or  surface  thermocouples.^ 

For  the  paint  tests,  the  wind  tunnel  model  is  electrodeposited  with  copper  to  a  thick¬ 
ness  of  .02  cm.  The  padnt  is  made  from  materials  of  calibrated  melting  points  suspended 
in  an  inert,  nonflammable  liquid.  The  paint  is  applied  to  the  model  and  allowed  to  dry  to 
a  dull,  opaque  film.  When  heat  is  applied,  the  film  liquefies  at  its  temperature  rating  to 
within  an  accuracy  of  l^C.^  For  the  thermocouple  tests  the  model  is  usually  made  of  steel 
and  is  hollow.  Thermocouples  axe  spot  welded  to  the  inner  surface  of  the  model  and  the 
lead  wires  are  passed  out  of  the  model  thru  a  base  seal.  The  air  temperature  inside  the 
model  is  usually  at  the  initial  tunnel  ambient  value,  which  acts  as  an  insulator  to  prevent 
heat  loss  firom  the  inner  surface.  The  thermocouples  are  calibrated  at  equilibrium  when 
immersed  in  water  baths  of  various  temperatures.  Paint  and  thermocouple  tests  yield  the 
time  variation  of  the  model  wall  temperature,  dT^fdt^  during  the  tunnel  nm.  Local  heat 
transfer  rate  is  obtained  firom  these  data  assuming  that  the  diffusion  time  through  the 
model  wall  is  very  small  competed  to  the  tunnel  run  time  (about  10  sec.)  and  using  the 
“thin  wall”  heat  balance  relation,^ 

.  _  .  dTy, 

—  PwCp^lw  ^ 

where  Pw,Cp^,  and  axe  the  model  wall  density,  heat  capeurity  and  thickness,  respectively. 

Local  heat  transfer  rates  have  been  me€isiu:ed  along  a  flat  plate  52  inches  (1.32  m)  in 
length  and  aligned  to  the  freestream  flow  at  a  Mach  and  Reynolds  number  of  10.53  and  25 
million  (based  on  total  length),  respectively.^®  The  plate  wall  temperature  was  conditioned 
to  4-times  the  freestream  value  before  the  test.  Flow  visualization  tests  indicated  that 
botmdary  layer  transition  &om  laminar  flow  started  at  .23  m  from  the  leading  edge  with 
full  turbulent  flow  developing  at  .5  m.  The  INTERACT  code  is  setup  for  both  laminar 
and  turbtilent  botmdary  layer  solutions  with  transition  specified  at  a  discrete  point.  A 
boundary  layer  transition  model  is  not  included  in  the  code.  Therefore,  for  the  flat  plate 
calculation,  transition  was  specified  at  .23  m  firom  the  leading  edge. 

Figure  10  shows  a  comparison  between  heat  transfer  rate  measurements  along  the 
plate  and  the  INTERACT  code.  Since  the  plate  has  a  sharp  leading  edge,  the  stagnation 
point  value  is  theoretically  infinite.  The  computed  laminar  heating  rates  slightly  xmderpre- 
dict  the  data.  The  maximum  heating  rate  after  transition  to  turbulent  flow  is  reproduced 

M.S.,  "An  Esperimtntnl  Inwtifntipn  of  TnrMent  Bpnninrf  Lnftri  at  Hi§k  Mnek  Nnmktr  nni  Refnelds  Nnm- 
ktr,"  Comtll  Aennnntienl  Lnkerntgrf  (Cnltpnn),  Bnffnlo  NY,  Report  No.  AB~S07t-A-l,  Noo.  197t  (tot  nUo  NASA  CR- 
Htt47). 
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by  the  code  at  the  specified  location.  The  computed  gradient  for  turbulent  heating  along 
the  remainder  of  the  plate  compares  well  with  the  data.  The  turbident  heating  level  b 
overpredicted  by  the  code  and  is  probably  a  strong  function  of  the  turbulence  model. 
Overall,  with  the  exception  of  flow  transition,  the  INTERACT  code  simulates  flat  plate 
aerodynamic  heating  within  acceptable  design  accuracy.  It  should  be  noted  that  the  hyper¬ 
sonic  fis^  plate  problem  is  not  trivial  since  it  involves  strong  shock-wave/botindary-layer 
interaction  along  the  boimdary  layer  edge,  particularly  near  the  plate  leading  edge,^  and 
is  a  good  test  of  the  inviscid-viscous  interaction  scheme. 

Figure  11  shows  a  compsoison  between  measured  and  computed  heat  transfer  rate 
along  a  15-degree  cone  with  a  20%  nose  blimtness.  These  data  were  gathered  during  actual 
high  altitude  (>  200,000  ft)  flight  tests  at  Mach  number  10.6  and  a  Reynolds  number  of  .38 
million  (based  on  total  length).^^  The  cone  wall  temperature  was  measured  at  a  constant 
value  1.76-times  that  of  the  freestream  static  value  before  data  acquisition.  The  data  have 
been  normalized  by  the  stagnation  point  heat  transfer.  The  calculation  overpredicts  heat 
transfer  along  the  blimt  nose  but  reproduces  the  constant  level  of  heating  along  the  cone 
section  of  the  body.  For  this  simulation,  modified- Newtonian  inviscid  theory  and  laminar 
boundary  layer  modeling  were  employed.  The  discontinuities  in  the  otherwise  smooth 
distribution  o'  heat  transfer  are  caused  by  rapid  boundary  layer  growth  near  the  stagnation 
point  and  body  slope  change  at  the  cone  jimction  {x/D  =  .06).  This  problem  is  also 
dominated  by  strong  shock- wave/botindary-layer  interaction.  Overall,  the  INTERACT 
code  predicts  the  heating  rate  to  within  design  accuracy. 

Figures  12-14  show  comparisons  between  computed  and  measured  Stanton  number 
(non-dimensional  heat  transfer,  see  Equation  66)  along  a  5-degree  blunt  cone,  cylinder 
midsection,  and  5-degree  flare  afterbody  projectile  as  determined  in  a  wind  tunnel  by 
RARDE.^  The  freestream  Mach  number  and  Reynolds  number  were  5.2  and  1.7  million 
(based  on  total  body  length),  respectively.  Both  zero  yaw  and  small  yaw  cases  were  tested. 
Some  tests  were  also  performed  for  Re  =  .65  and  3.1  million.  The  body  surface  temperature 
was  conditioned  at  .51-times  the  freestream  static  value  of  1482  R  before  the  tunnel  was 
started. 

Figure  12  shows  that  the  INTERACT  calculations  reproduce  the  zero-yaw  measured 
data  very  well.  Heating  levels  on  the  cone,  cylinder,  md  flare  sections  are  predicted. 
Figure  13  shows  a  comparison  between  the  computed  and  measured  zero-yaw  values  and 
the  transverse  plane  (^  =  90®)  values  at  1.25®  angle-of-attack  (a).  Since  the  code  does  not 
accovmt  for  cross  flow,  both  calc\ilated  distributions  are  identical  (see  Equations  56-57). 
The  measxired  data  show  a  cross  flow  effect  on  the  heating  rates  such  that  the  discrepancy 
between  measured  and  computed  rates  is  <  16%.  Cross  flow  effects  slightly  increase 
aerodynamic  heating  on  the  flare. 

Figure  14  shows  a  comparison  between  the  computed  and  measiired  heating  rates  in 
the  windward  (^  =  0®)  and  leeward  (^  =  180®)  planes  for  a  =  1.25®.  Once  again,  the  cross 
flow  effects  have  not  been  accoimted  for  in  the  INTERACT  code.  On  the  cone  section,  the 
measured  trend  of  higher  heat  transfer  on  the  windward  side  of  the  projectile  is  reproduced 
by  the  code,  however  the  magnitude  of  the  difference  is  underpredicted.  On  the  cylinder 

Goedriek,  W.D.,  Li,  C.P.,  Ho%*ton,  C.K.,  Chiu,  P.B.,  uni  Olmtio,  L.,  “NnmerieuI  Compntution*  of  Oriiter  Flowfitlit 
ond  Luminor  Heoiint  fUtei,"  AIAA  Jenmol  of  Spoeeernft  und  RoektU,  Vol.  14,  No.  S,  Mop  1977,  pp.  tS7-t84. 
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section,  the  same  trend  applies  to  the  measured  data.  The  calculation  predicts  a  reversal 
in  which  higher  heating  rates  occiir  on  the  leeward  side  at  the  end  of  the  cylinder.  On  the 
flare  afterbody,  both  measured  and  computed  results  show  that  heat  transfer  is  somewhat 
higher  on  the  windward  side  of  the  projectile.  Flare  heating  rates  are  reproduced  very  well 
by  the  INTERACT  code. 

Table  3  shows  a  comparison  between  the  computed  and  measured  heating  rate  at  a 
point  on  the  5-degree  flare  for  three  values  of  the  Reynolds  niimber.  The  heat  transfer  rate 
on  the  windward  side  for  a  =  1.25**  is  normalized  by  the  rate  for  zero  yaw.  The  error  in 
prediction  is  1%  to  1.5%  for  smaller  Re  and  12%  for  large  Re  values.  This  illustrates  that 
cross  flow  effects  on  heat  transfer  on  the  flsire  become  increasing  important  with  increasing 
Re. 


Table  3.  EFFECT  OF  REYNOLDS  NUMBER  ON  FLARE  HEAT  TRANSFER 


Comparison  made  at  X/D  =  17.6 


g  =  1.25°, ^  =  0°  normztlized  by  ,  a  =  0° 


Re(L) 

Measured 

Computed 

%  error 

.65  million 

1.249 

1.268 

1.5 

1.7  million 

1.092 

1.103 

1.0 

3.1  million 

1.723 

1.511 

12.3 

Figure  15  shows  a  comparison  between  the  computed  and  measured  heating  rates 
along  a  10-degree  flare  on  the  same  projectile  model.  The  zero- yaw  and  a  =  1.25®  windward 
side  are  compeured.  The  trend  zmd  overall  magnitude  of  the  data  are  predicted  by  the  code. 
Lower  heating  rates  mezisured  near  the  cylinder-flare  junction  are  not  predicted. 

Figure  16  shows  a  comparison  between  the  computed  and  measured  Stanton  niunber 
along  a  15-degree  cone,  cylinder  midsection,  and  10-degree  flare  afterbody  projectile  bs 
determined  in  a  wind  tunnel  test.^^  The  freestream  Mach  aind  Reynolds  numbers  were 
4.98  and  7.9  million  (bzised  on  toteil  body  length),  respectively.  The  model  was  at  zero 
yaw.  The  wall  temperature  was  cooled  to  .17-times  the  freestream  stagnation  value  of  710 
R  or  precisely  the  freestream  static  value  before  the  test.  The  heating  rates  on  the  flare 
afterbody  are  underpredicted  by  about  14%.  This  discrepancy  may  be  due  in  psurt  to  the 
steep  flaire  angle  where  rapid  boundary  layer  thickening  occmrs  on  the  flare. 


VI.  CONCLUSIONS 


An  engineering  design  code  utilizing  viscous-inviscid  interaction  has  been  developed 
to  obtain  aerodynamic  zmd  heat  transfer  predictions  for  axisymmetric  projectiles  at  high 
speed.  Separate  solutions  for  inviscid  and  viscous  flowfields  are  coupled  by  an  iterative 
interaction  model.  Inviscid  methods  include  linearized  perturbation,  second-order  shock- 
expansion,  and  modified-Newtonian  theories.  The  boundary  layer  equations  Me  solved 

^^Sckae/tr,  J.W.,  tni  Ferfuton,  H.,  “hvettigaticn  of  Separation  and  Aaaeeiated  Heat  Tranafer  and  Preaanrt  DiatrUntian 
an  Cone-Cplinder-Flart  Confitnraiians  at  Mach  Five.”  ARS  Journal.  Val.  3t,  Na.  S,  Map  l9St,  pp.  76t-770. 
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using  a  second-order,  finite-difference  scheme  with  specified  edge  velocity  (standard)  or 
displacement  thickness  (inverse).  The  procedure  has  been  formulated  for  supersonic  and 
hypersonic  Mach  numbers,  small  angle-of-attack,  and  Reynolds  ntunbers  for  laminar  and 
turbulent  flow.  The  capability  of  computing  local  reverse  flow  regions  is  included.  Ap¬ 
proximate  methods  for  stagnation  point,  laminar,  and  turbulent  heat  transfer  analysis  are 
also  included. 

Computations  of  projectile  surface  pressure  distribution  for  Mach  numbers  2  thru  5 
tmd  8  have  been  compared  to  wind  tunnel  measurements  on  several  cone-cylinder-flare 
configurations.  Computations  for  surface  heat  trtinsfer  coefficient  have  been  compared 
to  resiilts  obtained  from  wind  tunnel  measurements  on  flat  plate,  blunt  cone,  and  cone- 
cylinder-flare  models  at  Mach  numbers  5  and  10.  These  tests  show  that  the  code  is  able  to 
predict  the  surface  pressure  and  heat  transfer  to  a  reasonable  degree  of  accuracy.  Trends 
in  pressure  and  heat  transfer  with  Mach  number,  Reynolds  number,  flare  angle,  and  small 
angle-of-attack  have  been  reproduced.  Pressure  distributions  along  the  projectile  surface 
have  been  integrated  to  yield  force  and  moment  coefficients  and  center  of  pressure  predic¬ 
tions.  These  predictions  reproduce  wind  tunnel  data  to  a  reasonable  degree  of  accurawry. 
Predicted  aerodynamic  heat  transfer  distributions  along  the  projectile  surface  can  provide 
input  to  heat  conduction  codes  that  compute  material  stresses  and  deformations  due  to 
thermal  loads. 

The  INTERACT  code  yields  boundary-layer  and  surface  profiles  and  heat  transfer 
rates  in  about  one  minute  of  interactive  run  time  on  minicomputers.  Grid  generation  is 
not  required  and  input  parameters  are  read  from  a  namelist.  The  code  includes  an  extensive 
dictionary  of  input  parameter  definitions,  and  issues  warnings  when  these  parameters  are 
improperly  used  in  an  individual  or  global  sense.  The  code  has  been  written  using  standard 
FORTRAN  77  amd  can  be  run  on  personal  computers.  INTERACT  is  highly  suitable  for 
projectile  design  studies  aimed  at  maximizing  aerodynamic  performance  while  minimizing 
aerodynaunic  heating.  In  addition,  INTERACT  cam  be  used  to  test  and  improve  various 
algebraiic  turbulence  models. 
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Figure  1.  Linearized  Perturbation  Theory  -  Maximum  Body  Angle  vs.  Mach  Number 
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Figure  3.  Variation  of  (Cp)cone/(C',,)wedge  with  LOGio(Afoo^")- 
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Figure  4.  Wind  Tunnel  Data  (Cp)  vs.  Inviscid  Flow  (Re  =  oo)  and  INTERACT;  4®  FI 
Mach  =  2 
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Figure  5.  Comparison  of  Wind  Tunnel  Data  (Cp)  with  INTERACT  Results:  4®  Flare 
Mach  =  2,4  ’ 
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Figure  6.  Comparison  of  Wind  Tunnel  Data  (Cp)  with  INTERACT  Results:  6“  Flare, 
Mach  =  2,4,5 
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Figure  7.  Comparison  of  Wind  Tunnel  Data  (Cp)  with  INTERACT  Results;  8®  Flare, 
Mach  =  2,4,5 
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Figure  8.  Comparison  of  Wind  Tunnel  Data  (Cp)  with  INTERACT  Results;  20"  Flare, 
Mach  =  3,8.1 
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Figure  9.  Comparison  of  Wind  Tunnel  Data  (Cp)  with  INTERACT  Results-  6“  Flare 
Mach  =  2,  a  =  1®  ’ 
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Figure  10.  Comparison  of  Wind  Tunnel  Data  (q^)  with  INTERACT  Results;  Flat  Plate, 
Mach  =  10.53 
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Figure  12.  Comparison  of  Wind  Tunnel  Data  (g^)  with  INTERACT  Results;  5®  Flare 
Mach  =  5.2  ’ 
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Mach  =  5.2,  Re  (D)  =  84000,  Re  (L)  =  1710000 


Figure  13.  Compzurison  of  Wind  Tunnel  Data  {q^)  with  INTERACT  Results;  5®  Flare, 
Mach  =  5.2,  a  =  0®  and  1.25® 


51 


Stanton  Number  (St^)  *  1000 


RARDE  5“  Flare 

Mach  =  5.2,  Re  (D)  =  84000,  Re  (L)  =  1710000 


T  -  T 


0  2  4  6  8  10  12  14  16  18  20  22 


Axial  Distance  From  Nosetip  (cal) 


Figure  14.  Coinp2U'ison  of  Wind  Tunnel  Data  (gw)  with  INTERACT  Results;  5®  Flare, 
Mach  =  5.2,  a  =  1.25" 
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Figure  15.  Compau-ison  of  Wind  Tunnel  Data  (g„,)  with  INTERACT  Results;  10®  Flare, 
Mach  =  5.2,  a  =  0®  and  1.25® 
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Figure  16.  Comparison  of  Wind  Tunnel  Data  (g„)  with  INTERACT  Results;  10®  Flzure, 
Mach  =  4.98 
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1.  BRL  Report  Number  NR;-TR-3119 _ gf  jfgpgjt  90 _ 

2.  Due  Report  Received _ 

3.  Does  tfiis  rqxm  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

ftrr  which  the  report  will  be  used.) _ 


4.  ^tedficaDy,  how  is  the  report  being  used?  (Infonmaiion  source,  design  data,  procedure,  source 
of  id^  etc.) _ 


S.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  cdidencks  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  tedmical  content,  fonntt.  etc.) 


Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD  Organization 

ADDRESS  _ 

Address 


City.  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  suplc  or  tape  closed,  and  mail.) 
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Director  I 
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Aberdeen  Proving  Ground,  MD  210*  -3066 

OFFIQAL  BUSINESS 

BUSINESS  REPLY  MAIlI 

FIRST  CLASS  PERMIT  No  0001,  APG,  MD  | 
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Director 

U.S.  Army  Ballistic  Research  Laboratory 
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